Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1519 1744 (N35W61 L=345)  1074    113        2 (1474/1699-1577/1779)   5.2 W52-W74 N32-N36
      2 1618 1600 (N24W57 L=341)  1085    227        0 (1517/1549-1724/1634) -26.6 W48-W66 N18-N27
      3 1013 1436 (N26W07 L=292)   779    249        1 ( 917/1368-1128/1473)   7.6 W00-W15 N23-N27
      4  318 1389 (N36E47 L=237)   582     78        1 ( 272/1368- 366/1408) -14.0 E53-E41 N34-N38
      5  722 1376 (N28E13 L=272)   305     71        0 ( 696/1343- 759/1395)  14.0 E14-E10 N25-N28
      6  552 1055 (N11E28 L=256)   208     34        0 ( 533/1038- 576/1078) -33.7 E29-E26 N09-N12
      7 1643  896 (S17W37 L=322)  1410    297        0 (1587/ 773-1675/1041)  60.3 W32-W40 S23-S09
      8 1521  940 (S12W29 L=313)   270     65        0 (1495/ 903-1548/ 968)  36.9 W26-W30 S14-S11
      9 1356  834 (S16W16 L=301)   296     80        0 (1337/ 791-1380/ 876)  53.1 W14-W18 S17-S13
     10 1724  797 (S25W44 L=328)   381    107        0 (1695/ 743-1752/ 839)  50.2 W41-W46 S28-S22
     11  747  651 (S16E24 L=260)   316     89        0 ( 716/ 610- 785/ 693)  33.7 E26-E20 S17-S13
     12  251  494 (S17E77 L=207)   925     69        2 ( 226/ 465- 280/ 521) -53.1 E82-E73 S18-S14
  

    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