Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1150 1597 (N41W13 L=293)   306     73        0 (1117/1571-1180/1623) -45.0 W10-W15 N39-N43
      2 1021 1440 (N31W01 L=281)   310     45        2 ( 995/1414-1043/1468)  63.4 W00-W02 N28-N32
      3 1761 1279 (N16W56 L=336)   576    112        1 (1711/1242-1806/1307) -26.6 W52-W60 N13-N17
      4 1507 1197 (N13W33 L=312)   655    154        1 (1444/1147-1561/1250) -41.2 W28-W36 N09-N16
      5  667  855 (S05E23 L=257)   438    126        0 ( 600/ 820- 711/ 877) -15.9 E27-E20 S07-S03
      6  266  841 (S06E56 L=224)   669     99        2 ( 221/ 820- 316/ 863) -11.3 E61-E50 S07-S05
      7  712  599 (S22E23 L=257)   337    146        1 ( 681/ 549- 739/ 654)  81.9 E25-E20 S25-S18
      8  308  605 (S22E58 L=221)   542    107        1 ( 277/ 574- 346/ 643) -38.7 E61-E55 S23-S19
      9  559  609 (S21E34 L=246)   257     88        0 ( 540/ 578- 586/ 641) -63.4 E35-E32 S23-S18
     10  621  595 (S22E30 L=250)   263     75        0 ( 596/ 566- 641/ 629)  45.0 E31-E27 S24-S19
     11  881  461 (S32E13 L=267)   269     15        0 ( 859/ 443- 904/ 479)  26.6 E14-E10 S33-S30
     12 1307  418 (S38W20 L=300)   822    144        1 (1277/ 351-1348/ 476)  73.3 W18-W22 S43-S33
  

    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