Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  786 1758 (N48E36 L=293)   580    109        0 ( 735/1725- 836/1781)  50.2 E39-E31 N44-N50
      2  327 1401 (N16E57 L=271)  2359    268        7 ( 243/1285- 428/1503)  56.3 E62-E51 N08-N23
      3 1819 1435 (N36W71 L= 40)   736    121        0 (1797/1383-1844/1484) -90.0 W69-W74 N32-N39
      4 1467 1363 (N27W27 L=355)   941    184        0 (1447/1277-1492/1442)  79.7 W25-W29 N21-N32
      5 1271 1199 (N14W14 L=342)  1422    328        4 (1132/1149-1394/1249)   3.4 W04-W21 N10-N18
      6  356  939 (S12E45 L=283)   273     83        1 ( 330/ 908- 377/ 974) -36.9 E47-E42 S13-S09
      7 1003  698 (S19W02 L=331)   709    350        1 ( 934/ 564-1059/ 838) -64.8 E03-W07 S27-S10
      8  379  733 (S25E44 L=285)   305     79        1 ( 367/ 700- 390/ 774) -90.0 E45-E42 S26-S21
      9 1651  533 (S23W55 L= 24)  2219    372        2 (1544/ 436-1750/ 663)  60.6 W48-W60 S30-S13
     10 1032  309 (S48W12 L=340)  1693    318        5 ( 973/ 184-1089/ 441) -85.6 W07-W17 S62-S36
  

    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