Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1164 1570 (N31W09 L=  2)   627    148        1 (1091/1519-1231/1616) -33.7 W04-W13 N27-N34
      2  410 1363 (N18E42 L=311)   533    111        0 ( 390/1314- 428/1419) -66.8 E44-E39 N14-N21
      3 1837 1324 (N17W60 L= 53)   508    102        0 (1800/1309-1866/1336)  -8.1 W55-W63 N15-N17
      4  266 1242 (N11E53 L=299)   630    192        1 ( 216/1219- 319/1291) -26.6 E59-E47 N09-N14
      5  552 1190 (N07E30 L=323)   230     71        0 ( 528/1148- 572/1230)  68.2 E31-E28 N04-N09
      6  864 1091 (N00E10 L=343)   254     55        1 ( 822/1068- 890/1109) -14.0 E12-E08 N00-N01
      7  855  915 (S10E10 L=342)   511    167        0 ( 829/ 824- 891/ 998)  73.3 E12-E07 S15-S05
      8 1756  802 (S15W51 L= 44)   446    146        0 (1701/ 787-1804/ 819)  -5.7 W46-W56 S16-S14
      9  300  668 (S24E53 L=299)  1426    282        3 ( 255/ 591- 328/ 787)  53.1 E59-E48 S28-S16
     10  573  649 (S26E31 L=322)   581     83        1 ( 547/ 612- 605/ 702) -80.5 E33-E28 S28-S22
     11 1157  633 (S27W09 L=  2)   408    103        0 (1133/ 580-1181/ 682) -71.6 W07-W10 S30-S24
     12  510  486 (S37E40 L=313)   354     42        1 ( 485/ 465- 539/ 512) -26.6 E42-E38 S38-S35
  

    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