Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  764 1622 (N38E27 L=308)  5707    785       12 ( 466/1518-1106/1770)  -7.0 E80-W01 N32-N46
      2 1220 1428 (N28W10 L=345)   422    106        1 (1185/1384-1256/1472) -56.3 W07-W12 N25-N31
      3 1027 1429 (N27E03 L=332)   216     72        0 (1017/1392-1039/1467) -78.7 E04-E02 N24-N29
      4  271 1435 (N20E65 L=270)   614     58        1 ( 255/1401- 298/1464) -76.0 E68-E61 N18-N22
      5  915 1372 (N22E10 L=325)   554    161        2 ( 887/1300- 955/1452) -63.4 E12-E06 N17-N27
      6 1628 1376 (N28W42 L= 17)   275     57        0 (1599/1360-1656/1390) -11.3 W39-W44 N27-N29
      7  190 1229 (N07E66 L=269)  2955    190        4 ( 154/1146- 220/1328) -65.6 E72-E60 N01-N12
      8 1815 1235 (N20W59 L= 34)   614     80        1 (1790/1190-1831/1277) -90.0 W56-W61 N17-N22
      9 1778 1162 (N15W54 L= 29)   549     88        0 (1732/1142-1817/1188)  -6.3 W49-W58 N13-N16
     10 1787  959 (N03W54 L= 29)   601    105        0 (1746/ 915-1837/ 983)  24.0 W50-W60 N00-N04
     11 1392  936 (S01W23 L=359)   226     28        0 (1372/ 916-1414/ 957) -18.4 W22-W25 S02-N00
     12  370  814 (S17E44 L=291)   386    106        0 ( 329/ 789- 421/ 847) -14.0 E47-E39 S17-S14
     13 1850  827 (S05W64 L= 39)   518     58        1 (1818/ 802-1874/ 843)   8.1 W60-W68 S06-S03
     14 1145  704 (S17W10 L=345)  1202    401        0 ( 994/ 642-1337/ 799)  33.0 W01-W21 S22-S09
     15  759  523 (S33E15 L=320)  1511    312        6 ( 626/ 430- 885/ 576)  29.7 E26-E05 S41-S27
     16  536  398 (S45E38 L=297)   494     68        1 ( 508/ 378- 575/ 419)  -8.1 E41-E33 S46-S43
  

    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