Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  889 1854 (N61W04 L=157)   668     74        2 ( 856/1821- 927/1879)  26.6 W00-W10 N57-N64
      2  993 1539 (N32W05 L=158)   487    151        0 ( 937/1481-1049/1587) -58.0 W01-W07 N27-N36
      3 1756 1268 (N05W54 L=206)   341     34        0 (1735/1245-1776/1292)  26.6 W50-W56 N04-N06
      4  361 1243 (N21E43 L=110)  1127    176        2 ( 297/1194- 439/1285)   9.5 E48-E36 N18-N23
      5 1565  887 (S15W33 L=185)  1443    388        0 (1524/ 714-1597/1052)  81.5 W30-W35 S25-S05
      6 1851  991 (S12W59 L=212)   658     85        1 (1803/ 966-1886/1026) -11.3 W53-W63 S13-S10
      7 1895  898 (S18W66 L=218)  1956    216        3 (1853/ 822-1927/ 991) -84.8 W59-W71 S23-S12
      8 1843  866 (S20W59 L=211)   626     36        2 (1819/ 836-1866/ 890)  21.8 W55-W61 S21-S18
      9 1699  600 (S35W48 L=200)   682    113        0 (1653/ 562-1747/ 635)  24.0 W43-W52 S37-S33
     10  571  583 (S21E37 L=116)   246     64        0 ( 553/ 548- 596/ 612)  38.7 E38-E33 S23-S19
     11 1703  434 (S49W60 L=212)   568     58        0 (1674/ 412-1730/ 454) -12.5 W52-W70 S50-S46
  

    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