Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1334 1406 (N27W24 L=301)  1018    277        1 (1251/1337-1426/1511) -59.7 W18-W29 N21-N34
      2  476 1378 (N28E40 L=238)   309    119        0 ( 427/1352- 527/1398) -11.3 E44-E34 N26-N29
      3  452 1208 (N17E39 L=239)   296    101        1 ( 425/1153- 483/1250)  71.6 E40-E36 N13-N19
      4 1182  892 (S04W09 L=287)   238     63        1 (1155/ 864-1216/ 917)  53.1 W07-W11 S06-S02
      5 1286  832 (S09W16 L=293)   502     89        2 (1252/ 788-1313/ 878)  78.7 W13-W17 S11-S05
      6 1713  803 (S14W48 L=326)  1172    180        2 (1653/ 722-1772/ 877)  42.0 W43-W53 S18-S09
      7 1620  799 (S13W40 L=317)   471    180        0 (1585/ 725-1639/ 873)  77.5 W37-W41 S17-S08
      8 1798  690 (S22W61 L=338)  3329    277        7 (1718/ 586-1862/ 790) -53.7 W51-W70 S30-S15
      9  300  610 (S21E60 L=218)   855    130        1 ( 281/ 555- 320/ 676)  53.1 E62-E54 S24-S16
  

    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