QGIS » QGIS Application

This is not reprojection issue, this is interface issue, but still a bug

1. Open QGIS, add a layer, go to project properties, CRS
2. Turn on OTF, choose Google Mercator from recent -> Note that Proj string DIDN'T change
----If you click Apply, OK now - the reprojection won't turn on.
3. Choose any other CRS -> Proj string changed -> choose Google Mercator again -> Proj string changed
Apply, OK
Now back to normal.

correction, this applies for any projection, not only Google Mercator

when you first open properties/crs, selection of any proj from recent doesn't change proj sttring and consequently doesn't reproject.

The vector reprojection seems fixed, the raster one is not, I'm testing it but it seems we need to reopen this or #4583

confirmed:

trasformazione diretta of
(832.453740, 1353.648851)
PROJ.4: +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +to +proj=tmerc +lat_0=39.66825833333333 +lon_0=-8.133108333333334 +k=1 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
Error: latitude or longitude exceeded limits

Giovanni Manghi wrote:

See also #4583
I confirm that now the message that pop ups in the canvas is slightly different but the problem is still there.
I also still confirm that with certain data (see #4583) there is a freeze.

Can't reproduce the problem with master. Is master still affected?

Giovanni Manghi wrote:

I just updated under both Linux and Windows and while under Windows seems ok, under Linux I still get the

forward transform of
(309.320618, 934.085281)
PROJ.4: +proj=longlat +datum=WGS84 +no_defs +to +proj=tmerc +lat_0=39.66825833333333 +lon_0=-8.133108333333334 +k=1 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
Error: latitude or longitude exceeded limits

message. You can try the attached raster that is in the 3763 CRS. As already described the problem pop ups only when reprojecting to a geographic CRS, like wgs84.

I'm probably missing something. How do you produce that? On current master on Linux I get:

1. Clean project (defaults to wgs84, otfr off)
2. Insert dem_clipped
3. canvas zooms to dem_clipped in EPSG:3763
4. switch otfr on
5. 2 errors about empty reprojected (canvas) extent.
6. Zoom to layer
7. Content appears in EPSG:4326

Hi Jurgen,
see the attached screencast, this is what I'm getting in master.

ps, in Windows it work as expected, no annoying warning pop

Giovanni Manghi wrote:

Hi Jurgen,
see the attached screencast, this is what I'm getting in master.

ps, in Windows it work as expected, no annoying warning pop

which version of proj.4 are you running?

Jürgen Fischer wrote:

Giovanni Manghi wrote:

Hi Jurgen,
see the attached screencast, this is what I'm getting in master.

ps, in Windows it work as expected, no annoying warning pop

which version of proj.4 are you running?

Version: 4.7.0-2~oneiric2

Giovanni Manghi wrote:

ps, in Windows it work as expected, no annoying warning pop

so you're not running valuetool on windows?

Jürgen Fischer wrote:

Giovanni Manghi wrote:

ps, in Windows it work as expected, no annoying warning pop

so you're not running valuetool on windows? ;)

Giovanni Manghi wrote:

damn me... I didn't had the value tool active so I didn't think it was installed... sorry for the noise.

2564e6a79 adds support for the QgsCsException, so that python plugins can also catch that exception. Now the popup comes from python and show where the exception wasn't catched.

The following patch makes valuetool ignore that exception - no errors popping up anymore.

diff -u ~/.qgis/python/plugins/valuetool/valuewidget.py /tmp/valuewidget.py 
--- /home/fischer/.qgis/python/plugins/valuetool/valuewidget.py    2012-01-28 22:35:33.578309041 +0100
+++ /tmp/valuewidget.py    2012-01-28 22:34:38.694036879 +0100
@@ -151,7 +151,11 @@
             pos = mapPos
             if not mapCanvasSrs == layerSrs and self.iface.mapCanvas().hasCrsTransformEnabled():
               srsTransform = QgsCoordinateTransform(mapCanvasSrs, layerSrs)
-              pos = srsTransform.transform(mapPos)
+              try:
+                pos = srsTransform.transform(mapPos)
+              except QgsCsException, err:
+                # ignore transformation errors
+                continue
             isok,ident = layer.identify(pos)
             if not isok:
                 continue

Jürgen Fischer wrote:

2564e6a79 adds support for the QgsCsException, so that python plugins can also catch that exception. Now the popup comes from python and show where the exception wasn't catched.

c0cbe53a also sets the project's CRS to the layer's CRS, when it zooms to the first layer's extent...

c0cbe53a also sets the project's CRS to the layer's CRS, when it zooms to the first layer's extent...

I think the users will find this a huge usability improvement. Thanks Jurgen.

c0cbe53a also sets the project's CRS to the layer's CRS, when it zooms to the first layer's extent...

Hi Jurgen, just tested this feature:

it works as the project CRS really changes after the first "zoom to layer" (example project in wgs84 and layer in a projected CRS). At the same I see that the scalebar still shows degrees. If I enter the project properties and pick again the new CRS and then click "ok", then the scalebar goes ok.

it works as the project CRS really changes after the first "zoom to layer" (example project in wgs84 and layer in a projected CRS). At the same I see that the scalebar still shows degrees. If I enter the project properties and pick again the new CRS and then click "ok", then the scalebar goes ok.

Hi Jurgen, let me know if you want me to file a separate ticket about this issue.