Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1544 1343 (N21W38 L=350)  2254    510        5 (1391/1231-1670/1428) -37.9 W27-W46 N12-N27
      2 1845 1285 (N14W66 L= 18)   759    116        4 (1819/1241-1872/1327) -78.7 W62-W69 N11-N16
      3  869 1230 (N18E09 L=303)   277    134        0 ( 825/1165- 912/1281) -58.0 E11-E06 N13-N21
      4  285 1081 (N10E53 L=259)   511    165        3 ( 237/1064- 333/1113)   5.7 E57-E47 N08-N11
      5  567 1054 (N08E29 L=282)   270     55        0 ( 542/1034- 588/1077)  26.6 E31-E27 N06-N09
      6 1882  945 (S07W66 L= 18)  1109    237        0 (1834/ 903-1927/ 975) -21.0 W59-W74 S10-S05
      7  664  910 (S01E23 L=289)   222     53        1 ( 632/ 887- 691/ 929) -33.7 E25-E21 S02-N00
      8 1545  750 (S16W33 L=345)  1148    309        2 (1444/ 649-1650/ 898)  45.0 W26-W41 S21-S07
      9  173  829 (S06E69 L=243)   748    126        1 ( 158/ 799- 204/ 874) -45.0 E71-E64 S08-S03
     10 1633  558 (S30W44 L=355)  1541    249        3 (1538/ 522-1710/ 605)  12.1 W35-W51 S32-S27
     11  466  532 (S25E44 L=268)   773    122        2 ( 420/ 500- 530/ 574) -35.5 E47-E39 S27-S22
  

    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