Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  861 1798 (N40E33 L=329)   988    141        2 ( 800/1767- 932/1821)  -4.4 E39-E26 N37-N41
      2 1163 1703 (N39E06 L=355)  1744    283        3 (1037/1644-1242/1752) -13.2 E16-W01 N35-N43
      3  359 1608 (N19E64 L=297)  1167    128        1 ( 325/1548- 405/1653)  90.0 E67-E61 N14-N22
      4 1446 1611 (N39W17 L= 19)   279     70        0 (1420/1586-1479/1627)  45.0 W15-W20 N36-N41
      5  528 1421 (N08E39 L=322)   321     81        0 ( 492/1389- 553/1455)  63.4 E40-E37 N06-N10
      6 1106 1364 (N15E02 L=360)   449    181        0 (1004/1342-1178/1376)   6.3 E07-W02 N12-N16
      7 1554 1335 (N23W28 L= 29)   255     63        0 (1529/1312-1584/1363) -21.8 W25-W30 N21-N24
      8  204 1133 (S12E61 L=300)  1640    287        3 ( 123/1023- 277/1233)  39.5 E71-E54 S19-S05
      9  258 1087 (S15E54 L=307)   382     90        0 ( 237/1041- 280/1120)  59.0 E56-E52 S17-S12
     10  324  987 (S20E47 L=315)   465    153        0 ( 310/ 921- 333/1071)  83.7 E47-E45 S23-S14
     11 1549  868 (S05W35 L= 36)   248     82        1 (1536/ 827-1564/ 910) -78.7 W34-W36 S07-S02
     12  280  770 (S34E53 L=309)   678    129        1 ( 234/ 729- 317/ 830)  32.5 E59-E47 S36-S29
     13 1319  614 (S24W28 L= 29)   904    231        3 (1288/ 503-1362/ 719) -63.4 W25-W32 S31-S17
     14 1711  659 (S12W55 L= 56)   590    110        3 (1684/ 625-1747/ 715) -21.8 W50-W60 S13-S09
     15 1555  482 (S26W51 L= 53)  1910    323        5 (1445/ 396-1680/ 577)  63.4 W45-W58 S33-S17
     16 1601  388 (S30W64 L= 65)  1145    143        1 (1530/ 347-1664/ 430)  52.1 W59-W68 S34-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