| layer = QgsProject.instance().mapLayersByName('replacement_table')[0] | |
| test = QgsProject.instance().mapLayersByName('test')[0] | |
| #print (layer.name()) | |
| #print (test.name()) | |
| feat = layer.getFeatures() | |
| seat = test.getFeatures() |
The geometry type of a feature class can be specified in a .gfs file used by OGR to map from a GML document to a simple feature schema. The following numbers can be used to specify a specific geometry type:
-2147483647 Point25D
-2147483646 LineString25D
-2147483645 Polygon25D
-2147483644 MultiPoint25D
-2147483643 MultiLineString25D
-2147483642 MultiPolygon25D
| def vinc_pt(f, a, phi1, lembda1, alpha12, s ) : | |
| import math | |
| """ | |
| Returns the lat and long of projected point and reverse azimuth | |
| given a reference point and a distance and azimuth to project. | |
| lats, longs and azimuths are passed in decimal degrees | |
| Returns ( phi2, lambda2, alpha21 ) as a tuple | |
| Parameters: | |
| =========== | |
| f: flattening of the ellipsoid |
| from qgis.utils import iface | |
| from qgis.core import QGis | |
| from qgiscommand.command import command | |
| @command() | |
| def bbox(): | |
| layer = iface.activeLayer() | |
| def pbounds (Lextent,s): | |
| e = Lextent |
| SELECT airport_heliport.gid, airport_heliport.designator, ST_Buffer(st_makepoint(st_x(airport_heliport.geom), st_y(airport_heliport.geom)), 18520, 'quad_segs=32') AS geom FROM gis20141211.airport_heliport; |
| rec=0 | |
| def autoIncrement(): | |
| global rec | |
| pStart = 0 #adjust start value, if req'd | |
| pInterval = 50 #adjust interval value, if req'd | |
| if (rec == 0): | |
| rec = pStart | |
| else: | |
| rec = rec + pInterval | |
| return rec |
| [name]&vbnewline& round( [sqKM])&" Km<sup>2</sup>" |
| -- Calculate Km2 of a polygon stores in GCS WGS84 | |
| -- by default ST_Area returns the calculation in meters so the appropriate conversion was applied | |
| -- fir_area_20131205 was the example table used | |
| -- the spatial data was in a geom field so geom::geography was added | |
| SELECT ST_Area(geom::geography)/1000000 as sqKM from fir_area_20131205 |
| -- Test to compare the area calculation in degrees when data is stores as WGS84 in ESRI vrs POSTGIS | |
| select shape_area as ESRI_area, ST_Area (geom) as POSTGIS_area, (shape_area-(ST_Area (geom)))*100 as Percentage_Difference from fir_area_20131205 |
This is a non-technical reading list for technical people.
This is a list of software you should read like a novel.