Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  924 1822 (N45E28 L=163)   434     60        0 ( 890/1806- 956/1840)  26.6 E30-E25 N42-N46
      2  704 1713 (N31E38 L=153)   751    151        1 ( 634/1673- 757/1743)  -4.8 E43-E31 N29-N34
      3  601 1658 (N26E43 L=148)   560     85        1 ( 573/1614- 625/1702) -60.3 E44-E40 N22-N29
      4 1691 1559 (N43W44 L=235)  1129    126        2 (1650/1511-1744/1608) -29.1 W38-W49 N39-N45
      5  288 1530 (N14E66 L=125)  1981    115        5 ( 255/1473- 319/1584) -81.9 E73-E60 N10-N18
      6  507 1553 (N18E45 L=146)   300     73        1 ( 475/1532- 537/1581)  56.3 E47-E43 N15-N19
      7  347 1509 (N13E58 L=133)   552     72        3 ( 321/1471- 379/1539)  90.0 E60-E55 N10-N15
      8 1585 1368 (N25W31 L=221)   672    142        0 (1534/1297-1638/1404) -32.0 W26-W34 N21-N27
      9 1793 1266 (N24W50 L=241)   760    141        0 (1744/1203-1834/1313) -26.6 W44-W54 N20-N26
     10  174 1044 (S16E64 L=127)  4297    327       11 ( 117/ 920- 222/1191)  76.0 E72-E58 S23-S07
     11  252 1077 (S14E54 L=137)   763    116        1 ( 223/1039- 285/1120)  51.3 E57-E51 S16-S11
     12  317  968 (S20E47 L=144)   278     63        0 ( 293/ 951- 349/ 984)  21.8 E49-E44 S21-S18
     13  533  776 (S29E29 L=162)   975    282        3 ( 418/ 738- 669/ 811)   5.7 E38-E18 S31-S27
     14 1713  626 (S15W57 L=248)  1489    227        5 (1666/ 565-1773/ 692)  78.7 W52-W62 S19-S09
     15  215  606 (S44E72 L=119)   909     68        1 ( 185/ 578- 242/ 634)   7.4 E87-E64 S45-S41
     16  397  448 (S54E52 L=139)  1363    161        1 ( 330/ 395- 468/ 489)  11.6 E77-E38 S58-S50
     17 1548  433 (S30W54 L=245)  2917    369        3 (1388/ 390-1716/ 485)  30.3 W41-W67 S36-S22
  

    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