Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1559 1674 (N29W59 L=185)   610     68        1 (1530/1651-1588/1696) -36.9 W56-W61 N26-N30
      2 1608 1545 (N20W55 L=180)  1014    244        1 (1581/1456-1638/1635)  54.0 W50-W60 N14-N25
      3 1781 1137 (S07W53 L=178)  1278    312        2 (1700/ 979-1820/1226) -82.9 W47-W55 S16-N00
      4 1850 1180 (S07W62 L=187)   464     68        1 (1836/1155-1872/1209)  71.6 W60-W64 S08-S05
      5 1636  649 (S32W32 L=157)  1404    272        2 (1526/ 550-1731/ 736)  18.4 W21-W40 S35-S28
      6 1122  617 (S19E04 L=121)   234    109        0 (1068/ 603-1178/ 629) -33.7 E07-E01 S21-S17
      7  815  557 (S15E24 L=101)   254     63        3 ( 779/ 539- 848/ 579)   0.0 E26-E21 S16-S13
      8  925  472 (S23E20 L=105)   535    130        0 ( 867/ 442- 988/ 524) -63.4 E22-E16 S26-S18
      9  677  439 (S20E38 L= 87)   699    180        2 ( 607/ 365- 722/ 519) -85.6 E39-E35 S26-S13
     10  516  396 (S20E52 L= 73)  1812    384        5 ( 474/ 320- 582/ 454)  48.0 E59-E46 S25-S15
  

    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