Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  802 1532 (N33E10 L=108)   527    100        1 ( 761/1486- 846/1571) -39.8 E12-E06 N29-N36
      2  167 1302 (N25E62 L= 55)  1456    206        2 ( 147/1231- 186/1385) -84.3 E65-E59 N20-N30
      3 1136 1216 (N09W09 L=126)   260     79        0 (1091/1182-1179/1247) -45.0 W06-W10 N06-N11
      4 1833 1022 (S08W57 L=174)   369     79        0 (1811/ 983-1855/1059)  45.0 W54-W59 S10-S06
      5 1894  855 (S19W63 L=181)   758    155        1 (1851/ 820-1929/ 899)  26.6 W58-W68 S20-S16
      6  307  715 (S12E53 L= 65)   659     97        1 ( 267/ 690- 351/ 745) -23.2 E56-E49 S13-S10
      7 1098  504 (S34E01 L=117)   406    141        1 (1057/ 443-1122/ 571)  71.6 E04-W01 S38-S29
      8 1500  504 (S39W30 L=148)   690    141        2 (1441/ 457-1562/ 571)  37.9 W25-W34 S41-S34
      9  312  474 (S27E62 L= 55)  1878    277        3 ( 283/ 404- 340/ 536)  52.1 E68-E58 S31-S22
     10 1409  459 (S41W23 L=141)   274     66        0 (1392/ 440-1439/ 488) -45.0 W21-W25 S42-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