Bug report #8742

Reproject layer does not work with custom CRS

Added by Pedro Venâncio over 10 years ago. Updated almost 9 years ago.

Status:Closed
Priority:Normal
Assignee:-
Category:Processing/Core
Affected QGIS version:master Regression?:No
Operating System: Easy fix?:No
Pull Request or Patch supplied:No Resolution:fixed/implemented
Crashes QGIS or corrupts data:No Copied to github as #:17456

Description

Reproject layer works ok with QGIS default CRS, but doesn't do any reprojection with custom CRS.

I'm with QGIS 2.0.1 from OSGeo4W, doing transformations with grid files.

Some custom CRS I've tested:

+proj=tmerc +k=1 +lon_0=1 +no_defs +x_0=200000 +ellps=intl +nadgrids=C:\\proj_lib/ptLX_e89.gsb +units=m +lat_0=39.66666666666666 +wktext +y_0=300000 +pm=lisbon

+proj=tmerc +k=1 +lon_0=-8.131906111111112 +no_defs +x_0=200180.598 +ellps=intl +nadgrids=C:\\proj_lib/pt73_e89.gsb +units=m +lat_0=39.66666666666666 +wktext +y_0=299913.01

Cont_AAd_CAOP20121_Pinhel.zip (145 KB) Pedro Venâncio, 2013-10-02 04:58 AM

History

#1 Updated by Giovanni Manghi over 10 years ago

  • Assignee deleted (Victor Olaya)
  • Category changed from 94 to 44
  • Status changed from Open to Feedback
  • Target version changed from Version 2.0.0 to Future Release - High Priority

Hi Pedro,
is this a regression? What about on linux (it may use a different proj version)? Can you send me (if you can't attach them) sample data and grids?

#2 Updated by Pedro Venâncio over 10 years ago

Hi Giovanni,

On linux is the same.

I can't check if this is a regression, because on 1.8 I get another error when I try to change the CRS:

An error has occured while executing Python code:

Traceback (most recent call last):
  File "C:/Users/pedro/.qgis//python/plugins\\sextante\\gui\\CrsSelectionDialog.py", line 47, in okPressed
    self.authid = self.selector.selectedAuthId()
AttributeError: 'QgsProjectionSelector' object has no attribute 'selectedAuthId'

Python version:
2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)]

QGIS version:
1.8.0-Lisboa Lisboa, f440b60

Python path: ['C:/Users/pedro/.qgis//python/plugins\\\\sextante', 'C:\\\\Users\\\\pedro\\\\.qgis\\\\python\\\\plugins\\\\elevation', 'C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python', 'C:/Users/pedro/.qgis//python', 'C:/Users/pedro/.qgis//python/plugins', 'C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python/plugins', 'C:\\\\Program Files (x86)\\\\ArcGIS\\\\bin', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\bin\\\\python27.zip', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\Python27\\\\DLLs', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\Python27\\\\lib', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\Python27\\\\lib\\\\plat-win', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\Python27\\\\lib\\\\lib-tk', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\qgis\\\\bin', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\Python27', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\Python27\\\\lib\\\\site-packages', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\Python27\\\\lib\\\\site-packages\\\\PIL', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\Python27\\\\lib\\\\site-packages\\\\win32', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\Python27\\\\lib\\\\site-packages\\\\win32\\\\lib', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\Python27\\\\lib\\\\site-packages\\\\Pythonwin', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\Python27\\\\lib\\\\site-packages\\\\wx-2.8-msw-unicode', 'C:\\\\Users\\\\pedro\\\\.qgis\\\\python\\\\plugins\\\\mmqgis/forms', 'C:\\\\PROGRA~2\\\\Quantum GIS Lisboa\\\\apps\\\\qgis\\\\python\\\\plugins\\\\fTools\\\\tools']

I attach a sample shapefile. It's a selection of the portuguese caop, on EPSG:3763.

Thanks!

#3 Updated by Giovanni Manghi over 10 years ago

Pedro Venâncio wrote:

Hi Giovanni,

On linux is the same.

I can't check if this is a regression, because on 1.8 I get another error when I try to change the CRS:

so this is a SEXTANTE/processing only issue? the same tool in the "vector" menu works as expected?

#4 Updated by Giovanni Manghi over 10 years ago

  • Category changed from 44 to 94

I attach a sample shapefile. It's a selection of the portuguese caop, on EPSG:3763.

using your sample data, projections and grids the result in sextante/processing is a vector that has attributes but no geometries, correct? at least here it is what happens. The "save as..." option seems to work ok, do you confirm?

#5 Updated by Giovanni Manghi over 10 years ago

Giovanni Manghi wrote:

I attach a sample shapefile. It's a selection of the portuguese caop, on EPSG:3763.

using your sample data, projections and grids the result in sextante/processing is a vector that has attributes but no geometries, correct? at least here it is what happens. The "save as..." option seems to work ok, do you confirm?

I'm wrong, the output has geometries but it is given the wgs84 crs so it does not shows by default... it is needed to manually change its crs.

#6 Updated by Giovanni Manghi over 10 years ago

Giovanni Manghi wrote:

Giovanni Manghi wrote:

I attach a sample shapefile. It's a selection of the portuguese caop, on EPSG:3763.

using your sample data, projections and grids the result in sextante/processing is a vector that has attributes but no geometries, correct? at least here it is what happens. The "save as..." option seems to work ok, do you confirm?

I'm wrong, the output has geometries but it is given the wgs84 crs so it does not shows by default... it is needed to manually change its crs.

because the produced prj file is empty

#7 Updated by Pedro Venâncio over 10 years ago

using your sample data, projections and grids the result in sextante/processing is a vector that has attributes but no geometries, correct? at least here it is what happens. The "save as..." option seems to work ok, do you confirm?

Yes Giovanni, Save as.. works fine.

I'm wrong, the output has geometries but it is given the wgs84 crs so it does not shows by default... it is needed to manually change its crs.

because the produced prj file is empty

Here there are two problems:

a) Processing creates an empty .prj file when using custom CRS;

b) In addition to the blank .prj file, "Reproject layer" actually does not do any reprojection, the output geometry remain with same CRS of the input layer.

#8 Updated by Giovanni Manghi over 9 years ago

  • Status changed from Feedback to Open
  • Affected QGIS version changed from 2.0.1 to master

#9 Updated by Giovanni Manghi almost 9 years ago

  • Resolution set to fixed/implemented
  • Status changed from Open to Closed

seems to work fine in the latest master, please reopen if necessary.

#10 Updated by Giovanni Manghi almost 9 years ago

  • Category changed from 94 to Processing/Core

Also available in: Atom PDF