open-street-map-processing/generate_map_density_squares.pl at master · cpsolver/open-street-map-processing · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
#----------------------------------------------------------------
#     generate_map_density_squares.pl
#----------------------------------------------------------------

#  (c) Copyright 2023 by Richard Fobes at SolutionsCreative.com


#----------------------------------------------------------------
#  Initialization.

$unit_integer_per_degree = 10000000 ;
$latitude_integer_offset = ( 1080 ) * $unit_integer_per_degree ;
$longitude_integer_offset = ( -820 ) * $unit_integer_per_degree ;
$latitude_degrees_per_category = 0.333333 ;
$longitude_degrees_per_category = 0.333333 ;
$maximum_downward_vertical_category = int( ( 160 * $unit_integer_per_degree ) / $latitude_degrees_per_category ) ;
$maximum_rightward_horizontal_category = int( ( 360 * $unit_integer_per_degree ) / $longitude_degrees_per_category ) ;
$downward_pixels_per_category = 1 ;
$rightward_pixels_per_category = 1 ;
$downward_pixel_offset = 100 ;
$rightward_pixel_offset = 20 ;
$pixel_width = 2 ;
$pixel_height = 2 ;
$maximum_count = -1 ;
$maximum_latitude_integer = -999999999999 ;
$minimum_latitude_integer = 999999999999 ;
$maximum_longitude_integer = -999999999999 ;
$minimum_longitude_integer = 999999999999 ;


#----------------------------------------------------------------
#  Read each line of the input file.

while ( $input_line = <STDIN> )
{
    chomp( $input_line ) ;
#    print $input_line . "\n" ;


#----------------------------------------------------------------
#  Count the number of items in each category.

    if ( $input_line =~ /([0-9]{11}) +([0-9]{11})/ )
    {
        $latitude_integer = $1 ;
        $longitude_integer = $2 ;
        $downward_vertical_category = int( ( ( - $latitude_integer ) + $latitude_integer_offset ) / ( $unit_integer_per_degree * $latitude_degrees_per_category ) ) ;
        $rightward_horizontal_category = int( ( $longitude_integer + $longitude_integer_offset ) / ( $unit_integer_per_degree * $longitude_degrees_per_category ) ) ;
        $category = sprintf( "%d" , $downward_vertical_category ) . "_" . sprintf( "%d" , $rightward_horizontal_category ) ;
        if ( ( $downward_vertical_category > -1 ) && ( $downward_vertical_category <= $maximum_downward_vertical_category ) && ( $rightward_horizontal_category > -1 ) && ( $rightward_horizontal_category <= $maximum_rightward_horizontal_category ) )
        {
            $count_at_category{ $category } ++ ;
            if ( $count_at_category{ $category } > $maximum_count )
            {
                $maximum_count = $count_at_category{ $category }
            }
        }
#        print $latitude_integer . " " . $longitude_integer . " " . $category . "\n" ;


#----------------------------------------------------------------
#  For debugging.

        if ( $latitude_integer > $maximum_latitude_integer )
        {
            $maximum_latitude_integer = $latitude_integer ;
        }
        if ( $latitude_integer < $minimum_latitude_integer )
        {
            $minimum_latitude_integer = $latitude_integer ;
        }
        if ( $longitude_integer > $maximum_longitude_integer )
        {
            $maximum_longitude_integer = $longitude_integer ;
        }
        if ( $longitude_integer < $minimum_longitude_integer )
        {
            $minimum_longitude_integer = $longitude_integer ;
        }
    }


#----------------------------------------------------------------
#  Repeat the loop for the next input line.

}


#----------------------------------------------------------------
#  Write the results.

foreach $category ( keys( %count_at_category ) )
{
    $count = $count_at_category{ $category } ++ ;
    $opacity = sprintf( "%.2f" , ( 0.5 + ( 0.5 * ( log( $count ) / log( $maximum_count ) ) ) ) ) ;
    ( $downward_vertical_category , $rightward_horizontal_category ) = split( /_/ , $category ) ;
    $downward_vertical_position = ( $downward_vertical_category * $downward_pixels_per_category ) + $downward_pixel_offset ;
    $rightward_horizontal_position = ( $rightward_horizontal_category * $rightward_pixels_per_category ) + $rightward_pixel_offset ;
    $output_line = '<rect x="' . $rightward_horizontal_position . '" y="' . $downward_vertical_position . '" width="' . $pixel_width . '" height="' . $pixel_height . '" style="fill: #00b800; fill-opacity:' . $opacity . ';" />' ;
    print $output_line . "\n" ;
}


#----------------------------------------------------------------
#  For debugging, write high and low values.

#  result:  10782221712 09222273883 11793358189 08200015228
#  translation: +80 degrees and -80 degrees is latitude range, +/-180 degrees is longitude range

# print $maximum_latitude_integer . " " . $minimum_latitude_integer . " " . $maximum_longitude_integer . " " . $minimum_longitude_integer . "\n" ;