Форумы GIS-Lab.info

Код: Выделить всё

aPrj = aview.getprojection
distancedeg = Feature1.distance(Feature2)
distancemet = Feature1.returnprojected(aPrj).distance(Feature2.returnprojected(aPrj))

Код: Выделить всё

QgsFeature feature;
...
QgsDistanceArea da;
da.setSourceCrs( mQGisIface->mapCanvas()->mapRenderer()->destinationSrs().srsid() );
da.setProjectionsEnabled( true );

double length = da.measure( feature );

Код: Выделить всё

SEEK XCODEKART+substr(STXYFILE+space(12),1,12)+STR(STABSX,10,3)+STR(STABSY,10,3)
IF FOUND()
*** ЗАХВАТ КЛИЕНТА ВЫПОЛНЕН
    SELECT MAPY
    YYRECNO=KOD
    SELECT MAPX
    XXRECNO=KOD
    * ОПРЕДЕЛЕНИЕ КООРДИНАТ ОБЛАСТИ АНАЛИЗА
    QXBMPFILE=STXYFILE   
    TIK7=2
    DO WHILE TIK7<13
      X=SUBSTR(QXBMPFILE,TIK7,1)
      IF X="T"
         TTYBMPY=VAL(ALLTRIM(SUBSTR(QXBMPFILE,TIK7+1,12)))
         TTXBMPX=VAL(ALLTRIM(SUBSTR(QXBMPFILE,2,TIK7)))
         EXIT
      ENDIF
      TIK7=TIK7+1  
    ENDDO
    QXOTFILE="L"+ALLTRIM(STR(TTXBMPX-200,10,0))+"T"+ALLTRIM(STR(TTYBMPY-200,10,0))
    QXDOFILE="L"+ALLTRIM(STR(TTXBMPX+200,10,0))+"T"+ALLTRIM(STR(TTYBMPY+200,10,0))

    TIKX3=(TTXBMPX-XXRECNO)/100
    poutn3=(TTYBMPY-YYRECNO)/100
  
    XBMPX3=TIKX3*XMAHTAB+STABSX*XMAHTAB/100
    YBMPY3=poutn3*DELITY+STABSY*XMAHTAB/100
    KORXKLI=YBMPY3
    KORYKLI=XBMPX3
***  АНАЛИЗ ЧАСТИ КАРТЫ С КООРДИНАТАМИ  ОТ "+qXOTFILE+"   ДО "+qXDOFILE
     * ЦИКЛ ПО КВАДРАТАМ
     XOT1=TTXBMPX-200   
     YOT1=TTYBMPY-200
     XDO2=TTXBMPX+200
     YDO2=TTYBMPY+200
     * эти координаты абсолютные
     KORMINX=0
 	 KORMINY=0
	 MINDLINA=999999    
	 YYYTEK2=YOT1
	 DO whil YYYTEK2<YDO2+100
		XXXTEK2=XOT1
		DO whil XXXTEK2<XDO2+100
 	                    SKVADRAT="L"+ALLTRIM(STR(XXXTEK2,10,0))+"T"+ALLTRIM(STR(YYYTEK2,10,0))
	                    * ЦИКЛ ПО БАЗЕ STREETMA ПО ТЕКУЩЕМУ КВАДРАТУ
	                    SELECT STREETMA
	                    *INDEX ON CODEKARTA+BMPFILE1+STR(BMPX1,10,3)+STR(BMPY1,10,3) TAG STREETMA2
	                     SET ORDER TO  2
	                     SEEK  XCODEKART+SUBSTR(SKVADRAT+space(12),1,12)
	                     IF FOUND()          
  			                       TIKX=(XXXTEK2-XXRECNO)/100
			                       poutn=(YYYTEK2-YYRECNO)/100

    	                                * ЕСЛИ В КВАДРАТЕ НЕСКОЛОКО ОТРЕЗКОВ
	                            DO WHILE .NOT.EOF()
			                        IF CODEKARTA==XCODEKART.AND.BMPFILE1=SKVADRAT
                               * ПОИК МИНИМАЛЬНОГО РАССТОЯНИЯ ОТ ТОЧКИ ДО ОТРЕЗКА
                               * МЕТОД ДВИЖУЩЕЙСЯ ПРЯМОЙ
                                                       XBMPX1=BMPX1
 	                                       YBMPY1=BMPY1
	                                       XBMPX1=TIKX*XMAHTAB+XBMPX1*XMAHTAB/100
	                                       YBMPY1=poutn*DELITY+YBMPY1*XMAHTAB/100
	                                       XBMPX2=BMPX2
	                                       YBMPY2=BMPY2
	                                       XBMPX2=TIKX*XMAHTAB+XBMPX2*XMAHTAB/100
	                                       YBMPY2=poutn*DELITY+YBMPY2*XMAHTAB/100
 
                                       = SEEKLine(YBMPY1,XBMPX1,YBMPY2,XBMPX2,KORXKLI,KORYKLI)
	
     		                         SELECT STREETMA             
		   	   ENDIF
           		IF EOF()
	             EXIT
		        ELSE
	             SKIP
	            ENDIF
        	ENDDO
	      ENDIF
	     XXXTEK2=XXXTEK2+100    
	     IF XXXTEK2>XDO2
	       EXIT
	     ENDIF
             ENDDO  
	  YYYTEK2=YYYTEK2+100
	  IF YYYTEK2>YDO2
	     EXIT
	  ENDIF
    ENDDO


IF MINDLINA#999999
  * KORMINX=XTEK  
   * KORMINY=YTEK
   * ПРЕОБРАЗОВАТЬ ИЗ АБСОЛЮТЫХ КООРДИНАТ В КООРДИНАТЫ ЭКРАНА
...
...