Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1427 1798 (N54W62 L= 52)   756     53        2 (1404/1765-1454/1824)  14.0 W54-W71 N51-N55
      2 1339 1782 (N54W46 L= 36)   619     73        2 (1306/1750-1370/1806)   9.5 W39-W51 N51-N56
      3  457 1530 (N41E46 L=304)   804    121        2 ( 397/1505- 519/1551)   0.0 E54-E39 N38-N42
      4  826 1490 (N36E11 L=339)   511    153        0 ( 797/1420- 856/1550)  71.6 E13-E09 N31-N40
      5 1062 1454 (N32W06 L=356)  1066    234        3 ( 980/1414-1173/1517) -31.6 W00-W13 N28-N36
      6  570 1495 (N38E33 L=317)   282     53        1 ( 549/1465- 600/1513) -45.0 E35-E30 N35-N39
      7 1736 1360 (N19W57 L= 48)   369     68        1 (1720/1330-1757/1396)  68.2 W56-W59 N16-N21
      8  361 1350 (N28E50 L=301)   379     53        0 ( 341/1320- 387/1377)  45.0 E51-E47 N26-N30
      9 1087 1239 (N18W06 L=356)   655    227        1 (1027/1159-1171/1301)  35.0 W01-W11 N13-N21
     10 1628 1153 (N08W42 L= 32)   436     61        1 (1595/1115-1652/1184) -53.1 W39-W43 N05-N09
     11 1540 1095 (N05W34 L= 24)   536    146        1 (1460/1079-1599/1114) -11.3 W28-W38 N03-N06
     12 1782  788 (S17W55 L= 45)   479     65        1 (1759/ 768-1811/ 819)   9.5 W52-W58 S18-S14
     13 1156  558 (S26W06 L=356)   639    255        2 (1056/ 507-1266/ 635)  26.6 E01-W13 S28-S21
     14  565  507 (S26E38 L=313)   545    159        0 ( 503/ 446- 602/ 558) -63.4 E41-E35 S30-S22
     15  452  481 (S28E50 L=301)   827    211        0 ( 432/ 394- 500/ 557)  38.7 E56-E41 S34-S22
     16 1342  391 (S41W21 L= 11)   505    138        1 (1269/ 369-1403/ 412)   5.2 W14-W26 S41-S39
  

    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