Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1181 1420 (N25W21 L=  8)  1267    389        4 (1017/1376-1348/1487)  -9.9 W10-W33 N21-N28
      2 1600 1338 (N08W44 L= 31)   285     35        0 (1579/1323-1629/1357) -56.3 W42-W45 N06-N09
      3 1506 1199 (N03W32 L= 20)   664    134        1 (1462/1148-1545/1251)  20.6 W28-W36 N01-N05
      4  473 1125 (N26E32 L=315)  1450    273        2 ( 346/1088- 592/1159) -22.8 E43-E24 N21-N29
      5  185  956 (N21E63 L=285)  1610    219        2 (  95/ 927- 268/ 986)  -4.2 E80-E53 N18-N23
      6 1230  814 (S10W06 L=354)  1031    242        1 (1151/ 723-1294/ 918)  32.5 W00-W11 S13-S05
      7 1533  832 (S18W25 L= 12)   463    150        1 (1486/ 780-1593/ 864)  12.5 W20-W29 S19-S16
      8  383  808 (N09E46 L=302)   328     62        0 ( 363/ 772- 408/ 842) -68.2 E46-E44 N06-N11
      9  982  588 (S16E14 L=334)   557    152        0 ( 937/ 530-1036/ 643)  24.0 E18-E09 S19-S14
     10  718  567 (S11E30 L=317)   296    138        2 ( 692/ 520- 766/ 618) -90.0 E31-E28 S15-S07
     11 1718  574 (S42W39 L= 26)   701     84        1 (1682/ 547-1759/ 601)   0.0 W34-W43 S43-S39
     12 1582  461 (S45W24 L= 11)   354     62        1 (1564/ 436-1597/ 499)  63.4 W21-W25 S47-S42
     13  972  427 (S26E20 L=327)   378     73        1 ( 947/ 392- 994/ 467)  38.7 E22-E17 S28-S24
     14  494  350 (S21E59 L=288)  1158    131        1 ( 480/ 280- 513/ 418)  38.2 E67-E53 S26-S15
     15 1050  248 (S42E25 L=322)   704     89        0 (1033/ 203-1066/ 306)  56.3 E29-E21 S46-S37
     16 1331  228 (S56E06 L=341)  1470    214        4 (1243/ 198-1421/ 267) -15.5 E14-W04 S60-S50
     17  987  185 (S46E36 L=312)   637     80        2 ( 951/ 158-1021/ 207) -21.8 E40-E31 S48-S42
  

    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