Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  552 1580 (N42E38 L=326)   330     61        1 ( 530/1559- 574/1601) -14.0 E39-E35 N40-N43
      2  580 1359 (N26E30 L=333)  3695   1069       11 ( 285/1135- 866/1555)  -3.4 E62-E09 N11-N38
      3 1490 1510 (N33W38 L= 41)   348     90        0 (1449/1488-1534/1535) -31.0 W34-W40 N31-N34
      4  920 1454 (N32E06 L=357)   340    104        0 ( 872/1435- 980/1468) -14.0 E09-E01 N30-N32
      5 1316 1404 (N27W22 L= 25)   428    166        0 (1244/1372-1394/1424) -24.0 W16-W26 N24-N28
      6  254 1232 (N18E57 L=306)  1650    361        6 ( 223/1165- 292/1319) -45.0 E62-E52 N13-N22
      7 1540 1221 (N14W35 L= 39)   328     97        0 (1494/1189-1588/1259) -39.8 W32-W38 N11-N16
      8  225 1124 (N11E59 L=304)   522    109        1 ( 207/1068- 245/1167)  71.6 E61-E57 N07-N13
      9  400 1082 (N09E42 L=321)   269    116        0 ( 375/1024- 435/1126)  50.2 E44-E39 N05-N11
     10  710  650 (S19E22 L=342)   984    340        4 ( 596/ 541- 826/ 766)  36.4 E32-E13 S26-S11
     11 1407  742 (S15W24 L= 27)   261     62        2 (1371/ 727-1430/ 765)   0.0 W21-W25 S15-S13
     12 1767  636 (S25W56 L= 59)   443     82        2 (1742/ 622-1796/ 653)   0.0 W52-W59 S25-S23
     13  924  491 (S30E09 L=354)  1097    330        3 ( 822/ 382-1048/ 564)  35.8 E18-W00 S38-S25
  

    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