Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1003 1751 (N45E05 L= 99)   309     49        0 ( 978/1731-1028/1769)  26.6 E07-E02 N43-N46
      2 1160 1721 (N43W08 L=112)   390     62        0 (1127/1707-1200/1737)  -9.5 W05-W11 N41-N44
      3  569 1545 (N28E35 L= 69)   417     93        1 ( 535/1512- 610/1581) -35.5 E37-E30 N25-N30
      4  272 1527 (N27E66 L= 38)  1000     97        0 ( 253/1482- 289/1574) -25.0 E77-E60 N24-N31
      5 1166 1511 (N27W08 L=112)   217     71        0 (1138/1481-1196/1541) -36.9 W05-W09 N25-N29
      6  490 1488 (N24E40 L= 64)   325     77        0 ( 473/1446- 510/1522)  78.7 E40-E38 N21-N26
      7 1655 1429 (N24W44 L=148)   382     65        0 (1624/1399-1688/1451) -31.0 W41-W46 N22-N25
      8 1844 1367 (N22W64 L=168)  1279     88        1 (1821/1320-1877/1409)  45.0 W60-W71 N18-N25
      9 1347 1256 (N11W19 L=123)   197     30        1 (1329/1229-1360/1278) -71.6 W18-W20 N09-N12
     10 1606 1252 (N12W37 L=141)   260     45        1 (1573/1229-1628/1269)  11.3 W34-W39 N11-N13
     11 1713 1147 (N07W46 L=150)   370     64        3 (1689/1117-1732/1174) -26.6 W43-W47 N04-N08
     12  227 1136 (N02E57 L= 47)   488     93        0 ( 182/1124- 271/1153)   0.0 E62-E52 N01-N02
     13 1691 1016 (S01W44 L=148)  1093    223        0 (1654/ 963-1753/1091)  90.0 W41-W49 S04-N03
     14  185  951 (S09E62 L= 42)  1204    137        3 ( 163/ 876- 226/1011)  58.0 E65-E56 S13-S05
     15  891  974 (S07E08 L= 96)   200     66        0 ( 857/ 956- 926/ 989) -26.6 E09-E05 S08-S06
     16  779  788 (S19E15 L= 89)   472    114        1 ( 758/ 741- 799/ 850)  71.6 E16-E13 S21-S15
     17  338  821 (S17E47 L= 56)   294     32        1 ( 319/ 793- 355/ 845) -71.6 E49-E45 S18-S15
  

    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