Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1036 1822 (N54W13 L=348)   474     48        0 (1003/1803-1068/1838)  23.2 W09-W16 N52-N55
      2 1170 1804 (N50W25 L=360)   424     46        2 (1149/1780-1196/1825) -26.6 W22-W27 N48-N52
      3  606 1703 (N49E30 L=305)   593     46        2 ( 576/1672- 633/1733) -56.3 E33-E26 N45-N51
      4 1359 1561 (N29W29 L=  4)   333     84        1 (1337/1516-1386/1598)  39.8 W26-W32 N26-N31
      5  863 1496 (N29E06 L=329)   208     43        0 ( 836/1478- 890/1514)  14.0 E08-E04 N28-N30
      6  387 1447 (N31E44 L=291)   450     50        1 ( 358/1422- 416/1472)  21.8 E46-E41 N29-N32
      7  484 1407 (N27E34 L=301)   456     82        1 ( 443/1391- 533/1423)   0.0 E37-E30 N26-N28
      8  232 1174 (N15E56 L=279)   393     70        1 ( 209/1144- 254/1208)  36.9 E57-E53 N13-N17
      9  587  636 (S21E33 L=302)   333     91        3 ( 537/ 611- 623/ 659) -18.4 E36-E30 S22-S19
     10  412  608 (S21E48 L=287)   482     83        1 ( 390/ 573- 435/ 651) -90.0 E49-E45 S22-S17
     11  707  465 (S33E29 L=306)   691    192        1 ( 632/ 418- 787/ 518)  17.4 E37-E21 S35-S30
  

    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