Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1633 1490 (N21W55 L=135)   553     66        0 (1608/1456-1661/1522)  15.9 W52-W59 N18-N22
      2 1519 1436 (N21W42 L=123)   373     96        1 (1470/1412-1556/1469)   0.0 W37-W45 N19-N22
      3  196 1163 (N26E64 L= 17)  1036    125        4 ( 161/1112- 236/1218) -36.9 E71-E57 N23-N29
      4 1260 1195 (N12W17 L= 98)   278     81        1 (1219/1178-1297/1210) -26.6 W15-W19 N10-N13
      5 1827 1107 (S06W58 L=139)   877    240        2 (1760/1013-1869/1176) -56.3 W52-W62 S12-N00
      6 1376  928 (S05W19 L=100)   766    246        3 (1306/ 860-1435/1051) -71.6 W16-W22 S10-N02
      7 1709  975 (S10W44 L=125)   274     96        0 (1697/ 930-1720/1023)  78.7 W42-W45 S12-S07
      8 1761  923 (S15W48 L=129)   685     74        2 (1733/ 892-1797/ 961)  63.4 W45-W51 S16-S12
      9 1906  920 (S20W67 L=148)  1553    116        2 (1875/ 875-1944/ 955)   9.5 W62-W77 S22-S17
     10  501  872 (N06E36 L= 45)   256     97        1 ( 466/ 850- 550/ 904)   8.1 E38-E31 N04-N07
     11 1481  750 (S19W24 L=104)   325    135        1 (1469/ 678-1500/ 803)  69.4 W22-W25 S23-S15
     12  427  602 (S10E49 L= 32)   999    187        6 ( 367/ 533- 484/ 675) -63.4 E52-E45 S14-S04
     13  851  584 (S17E19 L= 62)   439    130        0 ( 830/ 521- 872/ 642)  82.9 E20-E17 S21-S13
     14  300  575 (S11E63 L= 17)   844     92        3 ( 277/ 533- 325/ 630)  63.4 E67-E60 S13-S07
     15 1563  571 (S33W29 L=110)  1034    150        2 (1490/ 541-1614/ 619)   0.0 W23-W33 S35-S30
     16  373  522 (S15E58 L= 23)   812     73        2 ( 347/ 481- 394/ 564)  63.4 E61-E55 S17-S11
     17  597  489 (S19E40 L= 41)   575    135        3 ( 572/ 419- 619/ 554)  48.4 E44-E36 S24-S15
     18  300  497 (S16E70 L= 11)  1823    114        3 ( 258/ 457- 342/ 543) -29.7 E77-E65 S19-S12
     19 1139  490 (S29E03 L= 78)   434     99        1 (1115/ 443-1168/ 536)  56.3 E05-W00 S31-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