Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  652 1775 (N58E24 L=153)   424     45        1 ( 620/1756- 678/1792) -12.5 E28-E19 N56-N60
      2 1749 1570 (N23W70 L=246)   648     79        0 (1721/1548-1779/1587) -26.6 W64-W75 N21-N24
      3  236 1484 (N40E67 L=110)   631     52        1 ( 205/1462- 269/1499)  -9.5 E77-E59 N38-N41
      4  670 1453 (N31E18 L=159)   252    103        0 ( 645/1418- 702/1498) -56.3 E19-E15 N28-N34
      5 1595 1289 (N08W40 L=217)   247     55        1 (1577/1264-1616/1318)  56.3 W38-W42 N06-N09
      6 1456 1227 (N06W29 L=206)   285    106        0 (1434/1175-1477/1277)  56.3 W26-W31 N03-N09
      7 1791  898 (S17W52 L=229)   401    161        0 (1761/ 832-1821/ 966)  53.1 W49-W55 S21-S13
      8 1210  786 (S17W08 L=185)   937    426        0 (1149/ 628-1256/ 959)  78.7 W04-W10 S26-S05
      9 1597  938 (S13W35 L=212)   225    123        0 (1527/ 926-1645/ 947)   0.0 W30-W39 S13-S11
     10 1662  880 (S17W40 L=217)   619     82        1 (1626/ 856-1711/ 910)   0.0 W37-W44 S18-S15
     11 1589  406 (S49W42 L=219)   831    126        1 (1543/ 382-1649/ 440) -10.6 W35-W51 S51-S45
  

    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