hi,
i tested the patch and it works ok (also on lat/long location). since nobody complained against this new feature, i committed
the patch to cvs
many thanks to Paul, the author of this patch
jachym
On Wed, Oct 18, 2006 at 10:00:53AM +0100, Paul Kelly wrote:
Oops, the line endings on that patch got corrupted - please look at this
one instead.
On Wed, 18 Oct 2006, Jachym Cepicky wrote:
>hallo,
>
>On Tue, Oct 17, 2006 at 10:56:58PM +0100, Paul Kelly wrote:
>>Hello Jachym
>>In order for the behaviour of r.out.gdal to be consistent with other
>>r.out.* modules, I think it would be more appropriate for it to first
>>use
>>r.resample to resample the map according to the current region (this
>>is
>>what the other modules do) into a temporary raster map and then export
>>that.
>
>could you point me to some example? how to create temporary raster file?
>i tryed something, see this patch
>
>On Wed, Oct 18, 2006 at 09:03:54AM +0200, Maciej Sieczka wrote:
>>Please, but only as an option. We can't remove any current behavior
>>unless it's a bug.
>>
>>In GRASS 7 region-sensitiveness in r.out.gdal should become the
>>default
>>however - for consistency with other r.out.* modules.
>>
>
>i added -r flag for "region sensitive"
>
>jachym
>
>
>--
>Jachym Cepicky
>e-mail: jachym.cepicky@centrum.cz
>URL: http://les-ejk.cz
>GPG: http://les-ejk.cz/gnupg_public_key/jachym_cepicky-gpg_public_key.asc
>-----------------------------------------
>OFFICE:
>Department of Geoinformation Technologies
>Zemedelska 3
>613 00, Brno
>Czech Republick
>e-mail: xcepicky@node.mendelu.cz
>URL: http://mapserver.mendelu.cz
>Tel.: +420 545 134 514
>
Index: r.out.gdal
RCS file: /grassrepository/grass6/scripts/r.out.gdal/r.out.gdal,v
retrieving revision 1.24
diff -u -r1.24 r.out.gdal
--- r.out.gdal 19 Aug 2006 12:52:24 -0000 1.24
+++ r.out.gdal 18 Oct 2006 08:07:59 -0000
@@ -22,6 +22,10 @@
#% key: l
#% description: list supported output formats
#%END
+#%flag
+#% key: r
+#% description: region sensitive output
+#%END
#%option
#% key: input
#% type: string
@@ -89,6 +93,7 @@
OUTPUT="$GIS_OPT_OUTPUT"
CREATEKEY="`echo "$GIS_OPT_CREATEOPT" | sed 's+,+ -co +g' | sed 's+^+-co +g'`"
METAKEY="`echo "$GIS_OPT_METAOPT" | sed 's+,+ -mo +g' | sed 's+^+-mo +g'`"
+REGION=""
if [ $GIS_FLAG_L -eq 1 ] ; then
gdal_translate | grep ':' | grep -v 'Usage' | grep -v 'The following'
@@ -99,6 +104,15 @@
#echo $INPUT $FORMAT $OUTPUT
#echo $GDAL_INPUT_FORMATS
+# region sensitive output
+if [ $GIS_FLAG_R -eq 1 ] ; then
+
+ TEMPRASTER=r.out.gdal_${INPUT}_$$
+ r.resample input=$INPUT output=$TEMPRASTER
+ INPUT=$TEMPRASTER
+fi
+
+
#fetch the input raster map
eval `g.findfile element=cell file=$INPUT`
if [ ! "$file" ] ; then
@@ -129,9 +143,14 @@
unset METAKEY
fi
+
#do it
CELLHD=`echo $file | sed 's+/cell/+/cellhd/+g'`
echo "Writing format: $FORMAT"
echo "Writing type: $TYPE"
gdal_translate -of $FORMAT -ot $TYPE $CREATEKEY $METAKEY $CELLHD $OUTPUT
+if [ $GIS_FLAG_R -eq 1 ] ; then
+ g.remove $TEMPRASTER
+fi
--
Jachym Cepicky
e-mail: jachym.cepicky@centrum.cz
URL: http://les-ejk.cz
GPG: http://les-ejk.cz/gnupg_public_key/jachym_cepicky-gpg_public_key.asc
-----------------------------------------
OFFICE:
Department of Geoinformation Technologies
Zemedelska 3
613 00, Brno
Czech Republick
e-mail: xcepicky@node.mendelu.cz
URL: http://mapserver.mendelu.cz
Tel.: +420 545 134 514
