Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1634 1562 (N19W58 L=310)  1229    195        1 (1567/1529-1715/1593) -30.3 W53-W65 N14-N22
      2 1878 1162 (S10W64 L=315)  1319    252        5 (1849/1109-1911/1207)  16.7 W59-W70 S13-S07
      3 1293 1026 (N01W15 L=267)   199     39        0 (1274/1000-1312/1058)  18.4 W13-W16 N00-N01
      4 1358  943 (S06W17 L=268)   369     80        1 (1308/ 922-1400/ 959) -31.0 W14-W19 S07-S04
      5 1612  773 (S23W30 L=282)  1140    247        1 (1504/ 712-1704/ 856)   6.7 W21-W38 S25-S21
      6 1827  775 (S32W50 L=301)   428     81        1 (1800/ 752-1849/ 798)  21.8 W46-W52 S32-S30
      7  789  451 (S21E30 L=221)   324    115        0 ( 760/ 420- 832/ 499) -80.5 E31-E28 S23-S17
      8  444  414 (S17E58 L=194)   684    196        4 ( 417/ 364- 475/ 458)  31.0 E62-E52 S20-S13
  

    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