[GRASS-user] r.proj -Segmentation fault

GRASS GRASS user
1

Hy!

I´m trying to re-project a GEOTIFF using r.proj. When executing my command,
I always get the message "Segmentation fault". I tried to find something in
the forum and I found the post:
http://osgeo-org.1803224.n2.nabble.com/GRASSLIST-4255-Segmentation-fault-while-using-v-proj-r-proj-tc1851719.html
Segmentation fault - GRASSLIST 4255 ... but I don´t really understand what
I have to do.

Before I was using v.proj and it worked without any problems.

My command:
r.proj in=tif_btr_84 out=dem30_tif location=reproj_tif_84
mapset=reproj_tif_84 method=bilinear resolution=30

Thanks in advance!

Cheers Tim

--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-proj-Segmentation-fault-tp6181331p6181331.html
Sent from the Grass - Users mailing list archive at Nabble.com.

2

On Thu, Mar 17, 2011 at 4:18 PM, TimNorwey <timmy_weyrer@live.at> wrote:

Hy!

I´m trying to re-project a GEOTIFF using r.proj. When executing my command,
I always get the message "Segmentation fault".

Which GRASS version do you use on what operating system?

I tried to find something in the forum and I found the post:
http://osgeo-org.1803224.n2.nabble.com/GRASSLIST-4255-Segmentation-fault-while-using-v-proj-r-proj-tc1851719.html

... the post is from 2002 :slight_smile: Perhaps 15 versions ago, so not much
related.

Segmentation fault - GRASSLIST 4255 ... but I don´t really understand what
I have to do.

Before I was using v.proj and it worked without any problems.

Note that v.proj and r.proj are somewhat unrelated.

My command:
r.proj in=tif_btr_84 out=dem30_tif location=reproj_tif_84
mapset=reproj_tif_84 method=bilinear resolution=30

From which to which projection do you want to reproject?

Are you sure that it is mathematically possible (certainly
r.proj should not just crash in case that it is not possible).

Markus

3

I´m using GRASS 6.4 on Mac OS X 10.6.6.

The GEOTIFF is in projection WGS84:
http://osgeo-org.1803224.n2.nabble.com/file/n6181437/WGS84_Ursprung.png

and I´m trying to re-project it to UTM 15:
http://osgeo-org.1803224.n2.nabble.com/file/n6181437/UTM15_Ziel.png

I don´t know if the re-projection is possible.

Thanks

Tim

--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-proj-Segmentation-fault-tp6181331p6181437.html
Sent from the Grass - Users mailing list archive at Nabble.com.

4

Hey.

Now I created a subset of the original GEOTIFF and tried it again, but it
doesn´t work. So I guess it has nothing to do with the size of my file.

The re-projection from WGS84 to UTM15 is generally possible. A friend of
mine tried it with some example data using ArcGIS, but I have no access to
ArcGIS.

Any ideas about other reason for that error?

Thanks in advance.

Tim

--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-proj-Segmentation-fault-tp6181331p6183278.html
Sent from the Grass - Users mailing list archive at Nabble.com.

5

On Fri, Mar 18, 2011 at 4:17 AM, TimNorwey <timmy_weyrer@live.at> wrote:

Hey.

Now I created a subset of the original GEOTIFF and tried it again, but it
doesn´t work. So I guess it has nothing to do with the size of my file.

Could you package the two locations with only this source map inside
for testing?

Markus

6

hey,

attached you can find the package.

I put the subset of my GEOTIFF in there. Did you meant that or did you mean
the original one?

Thanks,

Tim

http://osgeo-org.1803224.n2.nabble.com/file/n6184233/ReProjectWGS82toUTM15.zip
ReProjectWGS82toUTM15.zip

--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-proj-Segmentation-fault-tp6181331p6184233.html
Sent from the Grass - Users mailing list archive at Nabble.com.

7

On Fri, Mar 18, 2011 at 12:50 PM, TimNorwey <timmy_weyrer@live.at> wrote:

hey,

attached you can find the package.

I put the subset of my GEOTIFF in there. Did you meant that or did you mean
the original one?

I meant the locations.

Thanks,

Tim

http://osgeo-org.1803224.n2.nabble.com/file/n6184233/ReProjectWGS82toUTM15.zip

... File not found ....

ciao
Markus

8

OK ... the locations are in there too.

I can download the file, without any problems.

But I will attach it here again and hope that you can download it:
http://osgeo-org.1803224.n2.nabble.com/file/n6184337/ReProjectWGS84toUTM15.zip
ReProjectWGS84toUTM15.zip

Thanks

--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-proj-Segmentation-fault-tp6181331p6184337.html
Sent from the Grass - Users mailing list archive at Nabble.com.

9

On Fri, Mar 18, 2011 at 1:31 PM, TimNorwey <timmy_weyrer@live.at> wrote:

OK ... the locations are in there too.

Perfect, like this testing is easy.

I can download the file, without any problems.

But I will attach it here again and hope that you can download it:
http://osgeo-org.1803224.n2.nabble.com/file/n6184337/ReProjectWGS84toUTM15.zip

GRASS 6.4.1svn (utm15):~/grassdata > r.proj in=tif_test out=dem30_tif
location=reproj_tif_84 mapset=reproj_tif_84 method=bilinear
resolution=30
Input Projection Parameters: +proj=longlat +no_defs +a=6378137
+rf=298.257223563 +towgs84=0,0,0,0,0,0,0
Input Unit Factor: 1
Output Projection Parameters: +proj=utm +no_defs +zone=15 +a=6378137
+rf=298.257222101 +towgs84=0,0,0,0,0,0,0
Output Unit Factor: 1

Input:
Cols: 180 (180)
Rows: 180 (180)
North: 30.100000 (30.100000)
South: 30.050000 (30.050000)
West: -90.600000 (-90.600000)
East: -90.550000 (-90.550000)
EW-res: 0.000278
NS-res: 0.000278

Output:
Cols: 165 (166)
Rows: 188 (189)
North: 3332400.000000 (3332400.000000)
South: 3326760.000000 (3326730.000000)
West: 731280.000000 (731250.000000)
East: 736230.000000 (736230.000000)
EW-res: 30.000000
NS-res: 30.000000

Allocating memory and reading input map...
100%
Projecting...
100%
r.proj complete.

No problems here... the resulting map looks fine. My system:

uname -a
Linux north 2.6.33.7-desktop-2mnb #1 SMP Mon Sep 20 18:19:20 UTC 2010
x86_64 x86_64 x86_64 GNU/Linux

Someone on Mac needs to try.. Or Tim, you need to run it through
a debugger, see

http://grass.osgeo.org/wiki/GRASS_Debugging#Using_GDB

and post the full backtrace to a Trac ticket for analysis.

PS: I'll send the reprojected map to you.

Markus

10

No problems here... the resulting map looks fine. My system:

uname -a
Linux north 2.6.33.7-desktop-2mnb #1 SMP Mon Sep 20 18:19:20 UTC 2010
x86_64 x86_64 x86_64 GNU/Linux

Ok. Thank you very much for trying!

Someone on Mac needs to try.. Or Tim, you need to run it through
a debugger, see

http://grass.osgeo.org/wiki/GRASS_Debugging#Using_GDB

and post the full backtrace to a Trac ticket for analysis.

I´m currently working on a thesis to finish my Bachelor degree and I´m very
busy. As soon I have time, I will do that. I already took a look on it and
it seems to me a little bit complicated...

PS: I'll send the reprojected map to you.

Thank you I got it. It´s perfect. :slight_smile:

Cheers Tim

--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-proj-Segmentation-fault-tp6181331p6194907.html
Sent from the Grass - Users mailing list archive at Nabble.com.

11

No problems here... the resulting map looks fine. My system:

uname -a
Linux north 2.6.33.7-desktop-2mnb #1 SMP Mon Sep 20 18:19:20 UTC 2010
x86_64 x86_64 x86_64 GNU/Linux

Ok. Thank you very much for trying!

Someone on Mac needs to try.. Or Tim, you need to run it through
a debugger, see

http://grass.osgeo.org/wiki/GRASS_Debugging#Using_GDB

and post the full backtrace to a Trac ticket for analysis.

I´m currently working on a thesis to finish my Bachelor degree and I´m very
busy. As soon I have time, I will do that. I already took a look on it and
it seems to me a little bit complicated...

PS: I'll send the reprojected map to you.

Thank you I got it. It´s perfect. :slight_smile:

Cheers Tim

--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-proj-Segmentation-fault-tp6181331p6194915.html
Sent from the Grass - Users mailing list archive at Nabble.com.