Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1014 1591 (N34W04 L= 32)   486     98        1 ( 992/1540-1037/1643) -82.9 W02-W04 N30-N38
      2 1763 1427 (N19W57 L= 86)  1101    194        4 (1738/1369-1792/1491)  45.0 W53-W62 N15-N22
      3 1212 1212 (N08W13 L= 41)   206     48        0 (1175/1199-1240/1228)   0.0 W10-W14 N07-N09
      4  412 1047 (N04E39 L=349)  1101    249        1 ( 355/ 922- 451/1144) -72.9 E43-E36 S03-N10
      5  241 1078 (N07E54 L=334)   489    105        2 ( 210/1040- 272/1116)  45.0 E57-E51 N04-N09
      6  183  991 (N02E61 L=327)  2616    545       12 ( 121/ 922- 269/1084) -23.6 E70-E52 S01-N08
      7  270  772 (S11E55 L=333)   525     98        1 ( 248/ 734- 297/ 804)  38.7 E57-E52 S13-S09
      8  356  624 (S21E51 L=337)   429     56        0 ( 335/ 596- 377/ 656) -90.0 E52-E49 S22-S19
      9  955  398 (S42E11 L= 18)   821    170        1 ( 881/ 363-1039/ 456)  19.7 E17-E03 S44-S38
  

    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