Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  594 1610 (N52E28 L=252)   606    106        0 ( 542/1581- 651/1630) -12.1 E35-E21 N48-N53
      2  511 1602 (N52E39 L=241)   383     21        0 ( 494/1578- 531/1625) -76.0 E41-E35 N50-N54
      3 1678 1423 (N18W55 L=335)   887    217        0 (1595/1368-1752/1498) -50.7 W50-W61 N13-N24
      4 1004 1301 (N23W03 L=283)   365    133        2 ( 951/1268-1061/1342) -45.0 W00-W06 N20-N26
      5 1500 1247 (N11W34 L=314)   318    107        0 (1469/1223-1541/1289) -59.0 W32-W36 N09-N14
      6  870  808 (S04E13 L=268)   320    113        0 ( 815/ 788- 926/ 824)  -8.1 E15-E08 S05-S02
      7  271  725 (S03E59 L=221)   563    110        1 ( 251/ 665- 294/ 758)  80.5 E62-E55 S07-S01
      8 1021  451 (S30E10 L=270)  2430    494        5 ( 831/ 325-1166/ 624)  28.1 E28-W02 S37-S20
      9  735  480 (S24E29 L=251)   302     45        1 ( 714/ 452- 761/ 508) -90.0 E30-E27 S25-S21
     10 1520  461 (S40W29 L=309)   661     73        1 (1480/ 440-1558/ 486) -14.0 W24-W32 S42-S37
     11  502  383 (S28E54 L=226)   455     59        0 ( 488/ 347- 515/ 416)  56.3 E56-E51 S31-S25
  

    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