Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  889 1568 (N21E24 L=281)   475     98        1 ( 857/1523- 915/1611) -33.7 E27-E21 N17-N22
      2 1681 1302 (N30W34 L=340)  2119    446        2 (1641/1133-1754/1466) -54.2 W28-W42 N19-N40
      3 1576 1311 (N26W25 L=331)   266     92        0 (1556/1262-1601/1347)  90.0 W23-W26 N23-N29
      4 1873 1245 (N34W58 L=  3)   921    147        0 (1827/1176-1903/1299) -21.8 W50-W62 N30-N36
      5  179  994 (S29E63 L=242)  1735    201        2 (  87/ 970- 282/1018)   4.8 E86-E50 S30-S25
      6  845  833 (S22E05 L=301)   267     96        0 ( 827/ 790- 865/ 871) -68.2 E06-E02 S24-S19
      7 1717  841 (N04W48 L=353)   599     93        0 (1673/ 824-1771/ 861)  12.5 W44-W53 N02-N06
      8 1134  643 (S24W18 L=324)   614    240        0 (1022/ 618-1246/ 673)  36.9 W11-W23 S28-S19
      9 1473  549 (S19W43 L=349)   448     68        2 (1442/ 518-1495/ 583)  90.0 W41-W46 S20-S16
     10 1562  499 (S18W53 L=358)   329     53        1 (1535/ 479-1586/ 521)  45.0 W50-W54 S19-S16
     11 1551  420 (S22W58 L=  3)  1655    194        2 (1445/ 386-1616/ 466)  42.0 W51-W64 S27-S18
     12 1413  289 (S34W60 L=  6)   398     67        0 (1395/ 265-1438/ 320)  -5.7 W55-W65 S35-S32
  

    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