New answers tagged


Try this version (modified from gene's response): vl = QgsVectorLayer("Point?crs=", "temporary_points2", "memory") vl.dataProvider().addAttributes([ QgsField("X_coord", QVariant.Int), QgsField("Y_coord", QVariant.Int), QgsField("Azimuth", QVariant.Int) ]) new_features = [] for row in range(0,rows): for col in range(0,cols): rspnt = ...


1) If your pnt shapefile has only one element, use the iterator method next(): feat = pnt.getFeatures().next() pntxy = feat.geometry().asPoint() If the shapefile has many elements, use list comprehensions pntxy = [feat.geometry().asPoint() for feat in pnt.getFeatures()] # select pntxy[0] or pntxy[1] or... 2) instead of GDAL, why don't you use directly ...


The analysis is done by the Goddard Institute for Space Studies. You can find more information on the methodology at the GISS Surface Temperature Analysis (GISTEMP) website.


The correct date format (in this instance) is I think as advertised in your DescribeCoverage response like 2016-01-14T00:00:00.000‌​Z. To get at a date in your coverage you need to apply a 'Slice' along the Time Dimension. In WCS 2.0 this slicing is carried out by means of a domain subsetting operation. With GET/KVP requests the parameter name is SUBSET ...


Add cost evaluator, called Minutes Usage_Type Cost Units Minutes Build the route network Create a new route Use the Minutes evaluator as Impendance in the route properties

Top 50 recent answers are included