Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1056 1855 (N57E28 L= 35)  2153    217        5 ( 954/1825-1157/1880)  33.7 E36-E19 N50-N62
      2  845 1504 (N23E22 L= 41)   267     93        1 ( 815/1482- 891/1522)   0.0 E24-E18 N22-N24
      3 1490 1211 (N19W25 L= 89)   984    308        4 (1347/1169-1606/1252)   3.0 W15-W34 N16-N22
      4 1655 1101 (N16W39 L=103)   851    213        3 (1574/1025-1721/1149)  50.7 W34-W45 N09-N20
      5 1820 1044 (N16W56 L=120)   593     88        2 (1800/1003-1841/1084) -38.7 W53-W59 N14-N18
      6 1714  606 (S11W59 L=123)  1067    151        3 (1690/ 539-1740/ 673) -48.8 W55-W64 S15-S07
      7 1368  494 (S26W36 L=100)   937    189        1 (1295/ 466-1459/ 539)  41.6 W31-W41 S29-S20
      8 1071  269 (S50W29 L= 93)   361     12        0 (1053/ 248-1088/ 289) -21.8 W26-W31 S51-S47
      9 1267  216 (S48W52 L=116)  1208    144        1 (1213/ 188-1308/ 255)   4.1 W45-W61 S50-S44
  

    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