Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1400 1739 (N56W08 L=218)   868     85        1 (1350/1715-1437/1774)  24.0 W03-W12 N51-N60
      2 1170 1605 (N36E07 L=203)  1069    310        2 (1104/1463-1256/1714) -35.8 E15-W03 N28-N43
      3 1670 1495 (N44W38 L=248)   828     83        1 (1625/1461-1703/1522) -14.0 W33-W43 N41-N46
      4  325  933 (S24E46 L=164)   314     73        0 ( 310/ 901- 340/ 971)  90.0 E47-E44 S26-S21
      5 1495  779 (S05W35 L=244)   258     42        1 (1479/ 757-1521/ 801) -45.0 W32-W36 S06-S04
      6  382  552 (S49E42 L=167)   796     79        1 ( 351/ 515- 412/ 587)  45.0 E47-E37 S52-S45
      7  420  487 (S53E39 L=170)   468     78        0 ( 395/ 457- 456/ 510)  33.7 E43-E33 S56-S50
      8 1486  384 (S28W54 L=264)  2705    438        2 (1264/ 318-1638/ 450)  36.5 W40-W65 S37-S20
      9  511  403 (S58E29 L=181)   437     43        1 ( 485/ 381- 532/ 425)  45.0 E32-E25 S60-S55
  

    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