Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  455 1670 (N59E56 L= 30)   587     91        1 ( 415/1649- 496/1684)   8.7 E63-E47 N57-N60
      2 1590 1411 (N19W47 L=132)   752    238        0 (1536/1362-1661/1449) -45.0 W41-W51 N14-N21
      3 1315 1220 (N13W21 L=106)   250     78        1 (1273/1198-1352/1240)  11.3 W18-W23 N12-N14
      4  263 1058 (N17E55 L= 30)   536     50        1 ( 240/1025- 290/1083) -71.6 E57-E52 N15-N18
      5 1858 1000 (S11W61 L=146)   769     69        1 (1830/ 973-1891/1034)   0.0 W57-W66 S12-S08
      6 1029  567 (S22E07 L= 78)   907    257        0 ( 936/ 497-1107/ 658)  26.6 E14-W00 S25-S17
      7 1503  556 (S32W26 L=111)   405    115        0 (1458/ 506-1544/ 604)  39.8 W22-W29 S34-S29
      8  552  564 (S16E39 L= 46)   461     95        1 ( 524/ 519- 594/ 601)  39.8 E42-E35 S18-S13
      9 1090  465 (S30E05 L= 80)   880    306        0 ( 967/ 394-1209/ 532)  17.5 E15-W04 S33-S27
     10 1550  498 (S38W31 L=116)   400     59        0 (1513/ 480-1584/ 520)   0.0 W27-W34 S38-S36
     11  661  424 (S27E37 L= 48)   622     88        1 ( 637/ 386- 689/ 466)  74.1 E38-E34 S30-S23
     12 1337  440 (S37W13 L= 98)   294     67        0 (1305/ 418-1376/ 461)  21.8 W09-W15 S38-S35
  

    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