[GRASS-dev] making r.out.gdal region sensitive

[Jachym Cepicky]

hi
On Wed, Oct 18, 2006 at 09:56:22AM +0100, Paul Kelly wrote:
> [...]
> What do you think?
> 
> Paul

i think, that you are much better coder :-)

thanks for the hints, i'll have a look at this

jachym

> 
> 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 at 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 at 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

> _______________________________________________
> grass-dev mailing list
> grass-dev at grass.itc.it
> http://grass.itc.it/mailman/listinfo/grass-dev

-- 
Jachym Cepicky
e-mail: jachym.cepicky at 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 at node.mendelu.cz
URL:    http://mapserver.mendelu.cz
Tel.:   +420 545 134 514
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 191 bytes
Desc: Digital signature
Url : http://lists.osgeo.org/pipermail/grass-dev/attachments/20061018/5b4139bb/attachment.bin