Hi Michael,
Thank you for reply.
Seems I am running Saga 2.08 64 bit (see
http://webgis.uster.ch/temp/saga_about_tool.png) - it is the one
provided with the OSGeo4W installer and which integrates with QGIS
processing.
Should I try a newer version? Version 2.1? Would that work better with
big files?
Here is the dialogue with my reclassify settings:
http://webgis.uster.ch/temp/saga_reclass_settings.png
and here the error message (not very useful, I know):
http://webgis.uster.ch/temp/saga_about_tool_reclassify_error_message.png
I will probably try to install a newer SAGA version (hopefully I can do
this in parallel with OSGeo4W).
My task is not so super urgent. I just want to try to do all with FOSS
software.
Thanks a lot so far,
Andreas
Am 2014-11-18 18:48, schrieb Michael Treglia:
Andreas,
You're running 64bit SAGA, right? (if not sure, click the question
mark on the menu bar, and click "About" - should be towards the top).
I actually like SAGA because I can throw pretty giant files at it.
Just to test out and make sure, I just loaded in a 10GB GeoTIFF with
integer values, and did a single value reclass that took about 3
minutes. You likely need ample RAM and disk space, but given that you
do, it should be feasible. [Watching the system resources, you can see
the respective amount of memory get taken up as you work with data in
SAGA]
The 4GB limit seems like it could be a 32 bit issue. I was using the
latest version of SAGA (2.1.3), 64 bit, downloaded from:
http://sourceforge.net/projects/saga-gis/files/SAGA%20-%202.1/SAGA%202.1.3/
[2]
If you're really stuck and need to knock this task out of the way - if
you can get me the data I can do the reclass for you.
hope that helps,
mike
On Tue, Nov 18, 2014 at 10:04 AM, Andreas Neumann
<a.neum...@carto.net> wrote:
Hi Carlos,
I reprojected to a different swiss projection system and the
expression still does not work. I tested various projections,
including UTM-32N. Very strange ...
There clearly seems to be a bug in the raster calculator. I will
open a bug with high priority.
Thanks for testing and helping!
Andreas
Am 2014-11-18 15:36, schrieb Carlos Cerdán:
Hi Andreas:
It is solved: you have to reproject your TEMP layer. I've tried
with
UTM-32N zone and raster calculator works fine.
Why raster calculator (and SAGA) doesn't work with your original
CRS?
Well, I don't know, perhaps devs group have the answer.
It was a good challenge.
Carlos
2014-11-18 4:37 GMT-05:00 Andreas Neumann <a.neum...@carto.net>:
Hi Carlos,
Thank you for your help!
I tried with
(TEMP@1=249)*255 + (TEMP@1<249)*TEMP@1 + (TEMP@1>249)*TEMP@1
My layer is called "temp". However, the resulting file contains all
0.
I tried in both QGIS 2.4 and 2.6.
Is there something wrong with the input file? It was generated from
a PDF with ghostscript. Here is the metadata:
-----
GDAL provider
GTiff
GeoTIFF
Datensatzbeschreibung
//gis/gis/mapserverdata/uep/release_2014_11/temp.tif
TIFFTAG_DATETIME=2014:11:17 09:50:08
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
TIFFTAG_SOFTWARE=GPL Ghostscript 9.10
TIFFTAG_XRESOLUTION=600
TIFFTAG_YRESOLUTION=600
Kanal 1
STATISTICS_MAXIMUM=255
STATISTICS_MEAN=247.18714741724
STATISTICS_MINIMUM=0
STATISTICS_STDDEV=37.198102288306
Dimensionen
X: 71575 Y: 70867 Kanäle: 1
Leerwert
*Leerwert nicht gesetzt*
Datentyp
Byte - Acht Bit vorzeichenlose Ganzzahl
Pyramidenübersichten
Räumliches Bezugssystem des Layers
+proj=somerc +lat_0=46.95240555555556 +lon_0=7.439583333333333
+k_0=1 +x_0=600000 +y_0=200000 +ellps=bessel
+towgs84=674.4,15.1,405.3,0,0,0,0 +units=m +no_defs
Layerausdehnung (in urspünglicher Projektion des Layers)
0.0000000000000000,-70867.0000000000000000 :
71575.0000000000000000,0.0000000000000000
Kanal
Kanal 1
Kanal Nr
1
Keine Statistik
Noch keine Statistik gesammelt
------------
Thanks again if you have any idea what may be wrong with my file or
the above listed raster calc statement.
Andreas
Am 2014-11-17 18:38, schrieb Carlos Cerdán:
Hi Andreas:
I think QGIS raster calculator can make this job. Conditional
syntax
is a little bit... different. The sentence that can do it is:
(TIFF=249)*255 + (TIFF<249)*TIFF + (TIFF>249)*TIFF
Where TIFF is your raster layer. The trick is: if the sentence in
parentheses is true, it gets value 1, else zero; so we have to
multiply by 255 to replace the 249 values. Next we have to add the
other values.
Good luck,
Carlos Cerdán
PD. Please let me know if i'm wrong with this suggestion
2014-11-17 12:02 GMT-05:00 Andreas Neumann <a.neum...@carto.net>:
Hi,
I have a very simple task (at least I thought it is simple), but am
not able to do it properly.
I have a grayscale tiff value where I would like to replace the
value of 249 with 255 - all other values should stay as they are.
I tried with gdal "nearblack", but nothing changed - I started a
separate email-thread about that at the gdal list.
Then I tried with the QGIS raster calculator - but it does not
support conditional syntax.
I tried with Processing and GRASS - but it would stop at 8% - the
file is quite big.
Finally I came across SAGA "Reclassify with grid value". The
dialogue is very complicated and hard to understand. After a while
it would stop with an error message telling me that Saga may not be
properly configured - but it says that it is properly installed.
Quite confusing ...
I am trying on Win7 64bit.
If the file wouldn't be so big I would try this task in The Gimp -
where it is very simple to replace one pixel values with another -
using "Select by color" and "fill with a specific color".
Any idea how I could achieve my task with GDAL, QGIS, GRASS or
Saga?
Thanks,
Andreas
_______________________________________________
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user [1] [1] [1]
Links:
------
[1] http://lists.osgeo.org/mailman/listinfo/qgis-user [1] [1]
_______________________________________________
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user [1] [1]
_______________________________________________
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user [1] [1]
Links:
------
[1] http://lists.osgeo.org/mailman/listinfo/qgis-user [1]
_______________________________________________
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user [1]
_______________________________________________
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user [1]
Links:
------
[1] http://lists.osgeo.org/mailman/listinfo/qgis-user
[2]
http://sourceforge.net/projects/saga-gis/files/SAGA%20-%202.1/SAGA%202.1.3/
_______________________________________________
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user
