Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1301 1607 (N33W37 L=322)  1343    269        7 (1186/1552-1422/1648) -45.0 W30-W43 N26-N38
      2 1739 1427 (N10W61 L=346)   604    106        2 (1698/1393-1775/1461) -51.3 W58-W63 N07-N12
      3  747 1305 (N30E11 L=274)   364    104        1 ( 703/1283- 807/1322) -23.2 E13-E06 N27-N31
      4 1490 1248 (N08W33 L=318)   378     69        1 (1462/1215-1514/1274) -63.4 W31-W35 N06-N10
      5 1683 1216 (N01W46 L=331)   515    121        1 (1642/1197-1738/1249) -23.2 W43-W51 S01-N03
      6 1851 1090 (S12W59 L=344)   485     55        1 (1830/1063-1872/1117) -90.0 W56-W62 S13-S10
      7 1259  690 (S18W06 L=291)   205     86        0 (1237/ 649-1277/ 732)  45.0 W03-W07 S19-S15
      8  389  637 (S03E51 L=233)   299     58        0 ( 364/ 616- 418/ 655) -45.0 E53-E49 S04-S01
      9  304  552 (S07E65 L=220)   507     63        2 ( 275/ 530- 334/ 579) -45.0 E68-E61 S08-S05
     10 1088  477 (S27E10 L=275)   225     62        0 (1061/ 458-1121/ 494)   0.0 E12-E07 S28-S26
     11  368  465 (S14E63 L=221)   628     61        1 ( 346/ 431- 393/ 487)  45.0 E67-E59 S16-S12
     12  727  268 (S34E46 L=239)  2699    329        7 ( 634/ 178- 831/ 403) -86.6 E55-E39 S42-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