Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1120 1789 (N56E07 L=168)   526      3        0 (1096/1768-1145/1808)  53.1 E10-E04 N53-N58
      2 1035 1780 (N53E15 L=160)   419     66        0 (1013/1750-1061/1801) -78.7 E17-E13 N50-N55
      3 1233 1771 (N57W05 L=180)   507     68        0 (1202/1752-1268/1788)  23.2 W02-W09 N54-N58
      4 1304 1759 (N57W13 L=189)   424     60        1 (1282/1725-1327/1779) -59.0 W11-W16 N53-N59
      5  973 1721 (N46E18 L=157)   359     65        0 ( 956/1689- 991/1754) -51.3 E20-E15 N43-N49
      6 1410 1731 (N56W26 L=201)   532     49        1 (1381/1708-1439/1748)   0.0 W22-W30 N54-N58
      7  675 1583 (N31E35 L=140)   368     76        1 ( 655/1549- 696/1616) -59.0 E37-E32 N28-N33
      8 1331 1467 (N33W15 L=190)   272     66        0 (1300/1450-1368/1484)   0.0 W12-W18 N31-N34
      9  277 1443 (N16E66 L=109)  1045    110        2 ( 249/1396- 311/1483) -45.0 E70-E62 N12-N18
     10 1738 1345 (N30W51 L=226)   455     66        0 (1706/1320-1771/1364)  18.4 W47-W56 N28-N31
     11  617  989 (S07E25 L=150)   413     76        3 ( 589/ 950- 646/1022) -36.9 E27-E22 S09-S04
  

    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