source: subversion/applications/utils/ooc/georeference.pl @ 30494

Last change on this file since 30494 was 20278, checked in by richard, 9 years ago

harder better faster stronger

File size: 908 bytes
Line 
1#!/usr/bin/perl -w
2
3        # Convert raw TIFFs to properly squared-up GeoTIFFs
4
5        foreach $fn (@ARGV) {
6                $fn=~s/\..+$//;
7                print "Processing $fn\n";
8                $/="\n";
9#               $/="\r\n";
10
11                $cmd='';
12                $min_e=9999999; $max_e=0;
13                $min_n=9999999; $max_n=0;
14                open (INFILE, "$fn.csv") or die "Couldn't open $fn.csv: $!\n";
15                while (<INFILE>) {
16                        $inline=$_; chomp $inline;
17                        ($e,$a,$x,$n,$a,$y)=split(/,/,$inline);
18                        if ($e!~/\d+/) { next; }
19                        $cmd.="-gcp $x $y $e $n ";
20                        if ($e<$min_e) { $min_e=$e; }
21                        if ($n<$min_n) { $min_n=$n; }
22                        if ($e>$max_e) { $max_e=$e; }
23                        if ($n>$max_n) { $max_n=$n; }
24                }
25                close INFILE;
26
27                if (-e "${fn}_r.tiff") { unlink "${fn}_r.tiff"; }
28                system "gdal_translate $cmd -a_srs EPSG:27700 $fn.tif temp1.tiff";
29                system "gdalwarp -s_srs EPSG:27700 -t_srs EPSG:27700 -srcnodata '0 255 0' -dstnodata '0 255 0' -te $min_e $min_n $max_e $max_n temp1.tiff ${fn}_r.tiff";
30                unlink "temp1.tiff";
31        }
Note: See TracBrowser for help on using the repository browser.