Find_Lines_inSameGrid_by_LineID.m

function Lines_inSameGrid = Find_Lines_inSameGrid_by_LineID ( LineID )
    %% Find_Lines_inSameGrid_by_LineID function description:
    % Imput: LineID
    % Output:
    % Lines_inSameGrid: Line ID Set
    %
    % Example
    % Edge_ID_Set
    % Lines_inSameGrid = Find_Lines_inSameGrid_by_LineID ( 35261 )
    % Revision Notes:
    %       (10/04/14)
    % by shenghua chen
    % step 1: find the grid id set in which the line pass through
    % step 2: get all the line in the these grid id.
    
    %% step 1: find the grid id set in which the line pass through
    load Lines_in_GPS_Area_with_GridID.txt;
    kk = ismember(Lines_in_GPS_Area_with_GridID(:,1),LineID);
    TheRecord = Lines_in_GPS_Area_with_GridID(kk, :);
    GridID_num = TheRecord(2);
    GridIDSet_of_line = TheRecord(3:(GridID_num+2));
    
    
    %% step 2: get all the line in the these grid id.
    load GridID_LineID_near_GPS_Trajectory.txt;
    Line_in_grid_Set = [];
    for i_grid = 1:GridID_num
        % grid
        GridIDSet_of_line(i_grid);
        mark_gridID = ismember(GridID_LineID_near_GPS_Trajectory(:,1),GridIDSet_of_line(i_grid));
        GridID_lineID_record = GridID_LineID_near_GPS_Trajectory(mark_gridID,:);
        
        % some grid do not contain lines
        if size(GridID_lineID_record,1) == 0
            LineID_Set_inGrid_record = [];
        end
        
        % grid do contain lines
        if size(GridID_lineID_record,1) > 0
            LineNum = GridID_lineID_record(2);
            LineID_Set_inGrid_record = GridID_lineID_record(3:(2+LineNum));
        end
        
        Line_in_grid_Set = union(Line_in_grid_Set, LineID_Set_inGrid_record);
    end
    Lines_inSameGrid = Line_in_grid_Set;
    
    
    %     %% visualizaiton
    %     for i_line = 1:size(Lines_inSameGrid,2)
    %         Visualization_road_link_by_line_ID ( Lines_inSameGrid(i_line) );
    %     end
    
end