[Qgis-developer] SAGA kernel missing CRS?

[Paolo Cavallini]

Il 13/07/2015 08:40, Victor Olaya ha scritto:
> Ok, i tested with your CRS definition and it works fine
> 
> I create a shapefile and add that CRS definitino in a prj file. I open
> that file and the CRS is correctly identified as EPSG:27700 - OSGB
> 1936 / British National Grid.
> 
> Then I calculate kernel density, and the output sdat file has a prj
> file with exactly that same CRS definition (this is what the fix that
> i added does, just moving the CRS definition from  input to output).
> Can you check that it is not happening like that in your case, opening
> the prj file of the sdat file and comparing?
> 
> Computing contour lines with the density layer works as expected, also
> using the same CRS definition.

>> PROJCS["OSGB 1936 / British National Grid",GEOGCS["OSGB
>> 1936",DATUM["OSGB_1936",SPHEROID["Airy
>> 1830",6377563.396,299.3249646,AUTHORITY["EPSG","7001"]],TOWGS84[446.448,-125.157,542.06,0.15,0.247,0.842,-20.489],AUTHORITY["EPSG","6277"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4277"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",49],PARAMETER["central_meridian",-2],PARAMETER["scale_factor",0.9996012717],PARAMETER["false_easting",400000],PARAMETER["false_northing",-100000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","27700"]]

Hi Victor,

further tests:
* in several occasions the prj is apparently wrong, as you can see above
- it looks more of a qpj, and in any case it is not recognized by qgis;
I had the same also for a grass command (during the workshop, I could
not do more debug)
*
processing.runalg("gdalogr:contour","/tmp/processing/a06b812a38094eea8f21468a15731962/TARGETOUTGRID.sdat",0.001,"deaths","",None)
also returns a shp without prj (this seems a genuine bug).

Thanks.
-- 
Paolo Cavallini - www.faunalia.eu
QGIS & PostGIS courses: http://www.faunalia.eu/training.html