Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1002 1873 (N60W32 L=225)  1732    227        2 ( 931/1805-1106/1935)   9.0 W18-W56 N55-N62
      2 1215 1739 (N42W35 L=228)   282     46        0 (1193/1715-1238/1761)  18.4 W32-W38 N40-N42
      3 1323 1726 (N38W43 L=235)   386     59        0 (1286/1712-1357/1738) -21.8 W39-W45 N36-N39
      4  244 1157 (N25E53 L=139)   681    168        0 ( 229/1082- 263/1243) -74.7 E55-E51 N19-N30
      5 1661 1018 (S11W40 L=232)   304     92        0 (1630/ 986-1693/1059) -59.0 W37-W41 S13-S08
      6  815  720 (S11E18 L=174)   292    119        0 ( 785/ 663- 841/ 773) -90.0 E19-E17 S14-S07
      7 1877  722 (S35W65 L=258)  1246     85        3 (1845/ 689-1912/ 773) -40.6 W59-W76 S37-S31
      8 1186  713 (S19W03 L=196)   225     67        1 (1154/ 701-1228/ 733) -14.0 W01-W05 S20-S17
      9  239  658 (S04E65 L=128)   432     61        0 ( 225/ 629- 257/ 687)  33.7 E67-E61 S06-S02
     10 1061  621 (S22E06 L=186)   263     94        1 (1019/ 594-1105/ 645)   9.5 E09-E03 S23-S20
     11  954  449 (S31E18 L=175)   220     61        0 ( 938/ 421- 970/ 484)  45.0 E19-E15 S32-S28
     12  882  351 (S36E27 L=165)  1240    307        2 ( 736/ 333-1033/ 381) -21.8 E36-E16 S40-S30
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y