Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1424 1588 (N25W43 L=212)  4435    715        8 (1127/1531-1708/1624) -24.9 W23-W65 N14-N34
      2  968 1420 (N28W08 L=177)   468    115        1 ( 924/1374- 997/1467)  -9.8 W05-W09 N24-N32
      3 1916 1170 (S13W68 L=237)  3711    357        8 (1871/1077-1959/1292)  77.9 W63-W73 S19-S05
      4 1736  789 (S29W40 L=209)  2832    510        4 (1586/ 609-1800/1018)  38.4 W25-W49 S37-S16
      5  686  911 (N07E22 L=147)   294     99        1 ( 648/ 868- 731/ 934)   8.1 E25-E18 N04-N07
      6  839  833 (S01E15 L=154)   254     99        0 ( 799/ 791- 863/ 879) -71.6 E16-E14 S04-N02
      7 1601  785 (S24W29 L=198)   673    195        0 (1540/ 686-1654/ 876)  36.0 W23-W34 S28-S20
      8 1307  624 (S26W06 L=175)   336    102        0 (1276/ 565-1331/ 677)  45.0 W02-W08 S28-S22
      9  619  630 (S08E34 L=135)   252     76        2 ( 599/ 602- 636/ 650) -90.0 E35-E32 S09-S06
     10  942  529 (S21E18 L=151)   670    267        0 ( 861/ 448-1053/ 626)  19.4 E25-E08 S24-S18
     11  375  416 (S16E64 L=105)  1177    121        1 ( 336/ 379- 414/ 447) -45.0 E67-E60 S18-S13
     12  846  360 (S30E31 L=138)   405     98        0 ( 805/ 331- 888/ 388)   8.1 E35-E28 S31-S27
  

    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