Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1285 1670 (N35W37 L=133)   840    160        1 (1217/1645-1374/1692) -24.4 W30-W42 N31-N37
      2 1551 1591 (N23W52 L=148)   467     79        0 (1516/1566-1589/1613)   6.3 W47-W56 N22-N24
      3  558 1384 (N34E23 L= 73)  1497    292        2 ( 441/1305- 675/1434) -31.4 E33-E15 N26-N38
      4 1098 1333 (N19W11 L=108)  1273    418        1 ( 917/1271-1262/1434)   0.0 W00-W24 N16-N21
      5 1749  969 (S17W47 L=143)   484    109        0 (1719/ 927-1783/1020)  29.7 W43-W50 S18-S14
      6 1903  870 (S27W66 L=162)   626     77        0 (1866/ 835-1936/ 903)   4.8 W59-W72 S28-S25
      7 1702  797 (S26W41 L=137)   622     83        2 (1666/ 760-1731/ 840)  31.0 W37-W43 S28-S23
      8  854  725 (S11E16 L= 80)   512    190        1 ( 838/ 639- 881/ 817)  56.3 E19-E12 S15-S06
      9 1385  667 (S26W15 L=111)   355     59        0 (1365/ 628-1407/ 698)  53.1 W12-W16 S28-S24
  

    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