経度緯度からメッシュ番号を算出する

• 緯度経度からメッシュ番号を算出するサブルーチン
• パッケージ名を「lat_long_to_mesh」とする
• lat_long_to_mesh::mesh($Longitude, $Latitude, $MESH_CODE)と呼ばれることで
• 経度/緯度/メッシュサイズをで入れると、指定されたメッシュ番号に変換して返す
• メッシュサイズ引数(1=1km/2=500m/3=250m/4=125m/5=100m)

#!/usr/bin/perl
	use strict;
	use warnings;
#---------------------------------------------------------------------------------------------------------------------
	#緯度経度からメッシュ番号を算出するサブルーチン
	#パッケージ名を「lat_long_to_mesh」とする
	#lat_long_to_mesh::mesh($Longitude, $Latitude, $MESH_CODE)と呼ばれることで
	#経度/緯度/メッシュサイズをで入れると、指定されたメッシュ番号に変換して返す
	#メッシュサイズ引数(1=1km/2=500m/3=250m/4=125m/5=100m)
	
package lat_long_to_mesh;
	
sub mesh {
    	my($Longitude_out, $Latitude_out, $MESH_CODE) = @_;
    	#print $Longitude_out,"-",$Latitude_out,"-",$MESH_CODE,"\n";

					#-------------------------緯度経度→メッシュ番号算出----------------------------------------------
					# 1st mesh(latitude)
					#・緯度×60分÷40分＝ｐ 余りａ ｐ＝第１次地域区画を示す数字の上２けた　の計算
					#
					my $mesh_P = int($Latitude_out * 60 / 40);
					my $mesh_a = ($Latitude_out * 60) - ($mesh_P * 40) ;
					#print "一次メッシュ(P=$mesh_P a=$mesh_a)\n";
					# 2nd mesh(latitude)
					#・ ａ÷５分＝ｑ 余りｂ ｑ＝第２次地域区画を示す数字の上１けた　の計算
					
					my $mesh_Q = int($mesh_a / 5);
					my $mesh_b = $mesh_a - ($mesh_Q * 5) ;
					#print "二次メッシュ(Q=$mesh_Q b=$mesh_b)\n";
					
					# 3tn mesh(latitude)
					#・ｂ×60秒÷30秒＝ｒ 余りｘ ｒ＝第３次地域区画を示す数字の上１けた　の計算
					
					my $mesh_R = int($mesh_b * 60 / 30);
					my $mesh_x = ($mesh_b * 60) - ($mesh_R * 30) ;
					#print "三次メッシュ(R=$mesh_R x=$mesh_x)\n";
					
					# 4tn mesh(latitude)
					#・ｘ÷15秒＝Ｒ２ 余りｘ２ Ｒ２＝第４次地域区画を示す数字の上１けた　の計算
					
					my $mesh_R2 = int($mesh_x / 15);
					my $mesh_x2 = $mesh_x - ($mesh_R2 * 15) ;
					#print "四次メッシュ(R2=$mesh_R2 x2=$mesh_x2)\n";
					
					# 5th mesh(latitude)
					#・ｘ２÷7.5秒＝Ｒ３ 余りｘ３ Ｒ３＝第５次地域区画を示す数字の上１けた　の計算
					
					my $mesh_R3 = int($mesh_x2 / 7.5);
					my $mesh_x3 = $mesh_x2 - ($mesh_R3 * 7.5) ;
					#print "五次メッシュ(R3=$mesh_R3 x3=$mesh_x3)\n";
					
					# 6rh mesh(latitude)
					#・ ｘ２÷3.75秒＝Ｒ３ 余りｘ３ Ｒ３＝第５次地域区画を示す数字の上１けた　の計算
					
					my $mesh_R4 = int($mesh_x3 / 3.75);
					my $mesh_x4 = $mesh_x3 - ($mesh_R4 * 3.75) ;
					#print "六次メッシュ(R4=$mesh_R4 x4=$mesh_x4)\n";
					
					# 100m mesh(latitude)
					# 緯度差3秒
					#・ｘ÷3秒＝Ｒ５ 余りｘ４ Ｒ４＝100mメッシュを示す数字の上１けた　の計算
					
					my $mesh_R5 = int($mesh_x / 3);
					my $mesh_x5 = $mesh_x - ($mesh_R2 * 3) ;
					#print "100mメッシュ(R5=$mesh_R2 x5=$mesh_x2)\n";
					
					# 1st mesh(longitude)
					#・経度－100度＝ｓ 余りｃ ｓ＝第１次地域区画を示す数字の下２けた　の計算
					#
					my $mesh_S = int($Longitude_out - 100);
					my $mesh_c = $Longitude_out -100 - $mesh_S ;
					#print "一次メッシュ(S=$mesh_S c=$mesh_c)\n";
					
					
					# 2nd mesh(longitude)
					#・ｃ×60分÷７分30秒＝ｔ 余りｄ ｔ＝第２次地域区画を示す数字の下１けた　の計算
					#
					my $mesh_T = int($mesh_c * 60 / 7.5);
					my $mesh_d = ($mesh_c * 60) - ($mesh_T * 7.5) ;
					#print "二次メッシュ(T=$mesh_T d=$mesh_d)\n";
					
					# 3th mesh(longitude)
					#・ｄ×60秒÷45秒＝ｕ 余りｙ ｕ＝第３次地域区画を示す数字の下１けた　の計算
					#
					my $mesh_U = int($mesh_d * 60 / 45);
					my $mesh_y = ($mesh_d * 60) - ($mesh_U * 45) ;
					#print "三次メッシュ(U=$mesh_U y=$mesh_y)\n";
					
					# 4th mesh(longitude)
					# 経度差は22.5秒
					#・ｙ÷22.5秒＝ｕ２ 余りｙ２ ｕ２＝第４次地域区画を示す数字の下１けた　の計算
					#
					my $mesh_u2 = int($mesh_y / 22.5);
					my $mesh_y2 = $mesh_y - ($mesh_u2 * 22.5) ;
					#print "四次メッシュ(u2=$mesh_u2 y2=$mesh_y2)\n";
					
					# 5th mesh(longitude)
					#・ｙ２÷11.25秒＝ｕ３ 余りｙ３ ｕ３＝第５次地域区画を示す数字の下１けた　の計算
					#
					my $mesh_u3 = int($mesh_y2 / 11.25);
					my $mesh_y3 = $mesh_y2 - ($mesh_u3 * 11.25) ;
					#print "五次メッシュ(u3=$mesh_u3 y3=$mesh_y3)\n";
					
					# 6th mesh(longitude)
					#・ｙ３÷5.625秒＝ｕ４ 余りｙ４ ｕ４＝第６次地域区画を示す数字の下１けた　の計算
					#
					my $mesh_u4 = int($mesh_y3 / 5.625);
					my $mesh_y4 = $mesh_y3 - ($mesh_u4 * 5.625) ;
					#print "六次メッシュ(u4=$mesh_u4 y4=$mesh_y4)\n";
					
					# 100m mesh(longitude)
					#・ｙ÷4.5秒＝ｕ５ 余りｙ５ ｕ５＝100mメッシュを示す数字の下１けた　の計算
					#経度差は4.5秒
					my $mesh_u5 = int($mesh_y / 4.5);
					my $mesh_y5 = $mesh_y - ($mesh_u2 * 4.5) ;
					#print "100mメッシュ(u5=$mesh_u2 y5=$mesh_y2)\n";
					
					#-------------組み立て-----------
					my $mesh1 = $mesh_P * 100 + $mesh_S;
					my $mesh2 = $mesh_Q *  10 + $mesh_T;
					my $mesh3 = $mesh_R *  10 + $mesh_U;
					my $mesh4 = oct('0b'.$mesh_R2.$mesh_u2) + 1;
					my $mesh5 = oct('0b'.$mesh_R3.$mesh_u3) + 1;
					my $mesh6 = oct('0b'.$mesh_R4.$mesh_u4) + 1;
					
					my $mesh_code_out = 0;
					
					if ($MESH_CODE == 1){$mesh_code_out = sprintf("%04d%02d%02d", $mesh1, $mesh2, $mesh3);} #1kmメッシュ計算
					elsif ($MESH_CODE == 2){$mesh_code_out = sprintf("%04d%02d%02d%01d", $mesh1, $mesh2, $mesh3, $mesh4);}#500mメッシュ計算
					elsif ($MESH_CODE == 3){$mesh_code_out = sprintf("%04d%02d%02d%01d%01d", $mesh1, $mesh2, $mesh3, $mesh4, $mesh5);}#250mメッシュ計算
					elsif ($MESH_CODE == 4){$mesh_code_out = sprintf("%04d%02d%02d%01d%01d%01d", $mesh1, $mesh2, $mesh3, $mesh4, $mesh5, $mesh6);}#125mメッシュ計算
					elsif ($MESH_CODE == 5){$mesh_code_out = sprintf("%04d%02d%02d%01d%01d", $mesh1, $mesh2, $mesh3, $mesh_R5, $mesh_u5);} #100mメッシュ計算
    	return $mesh_code_out;
}

---