Highways = shaperead('HE_HighwayBoundary_v3_3_1');

prof = geotiffinfo('GBOverview');

HighwaysLat = cell(length(Highways),1);
HighwaysLon = cell(length(Highways),1);

for i = 1 : length(Highways)
    highwaysX = Highways(i).X;
    highwaysY = Highways(i).Y;
    
    [HighwaysLat_temp, HighwaysLon_temp] = projinv(prof, highwaysX, highwaysY);
    
    Highways(i).X = HighwaysLat_temp;
    Highways(i).Y = HighwaysLon_temp;
    
    Highways(i).BoundingBox(1) = min(Highways(i).X);
    Highways(i).BoundingBox(2) = min(Highways(i).Y);
    Highways(i).BoundingBox(3) = max(Highways(i).X);
    Highways(i).BoundingBox(4) = max(Highways(i).Y);
end

save('HighwaysLatLon', 'Highways')