Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  787 1688 (N33E28 L=308)   662    202        0 ( 709/1616- 853/1737) -26.6 E36-E20 N28-N36
      2  592 1660 (N29E42 L=294)   397     98        1 ( 575/1620- 612/1709) -54.5 E44-E38 N26-N33
      3  666 1465 (N17E29 L=307)   316     90        1 ( 641/1414- 683/1498) -59.0 E31-E26 N13-N18
      4 1275 1327 (N15W11 L=347)   368    106        0 (1236/1291-1321/1365) -26.6 W08-W14 N13-N17
      5 1636 1293 (N19W37 L= 13)   384     59        2 (1601/1269-1663/1316) -11.3 W33-W38 N18-N21
      6  606  905 (S18E25 L=311)   465    122        0 ( 571/ 839- 646/ 951)  66.8 E27-E22 S21-S14
      7  207  850 (S24E60 L=275)   476     71        0 ( 177/ 819- 232/ 886) -38.7 E63-E57 S25-S21
      8  843  759 (S24E08 L=328)   576    189        1 ( 788/ 681- 879/ 851)  74.7 E10-E05 S29-S18
      9 1225  652 (S25W18 L=354)   195     53        0 (1196/ 633-1253/ 671)  45.0 W16-W19 S26-S23
     10 1712  644 (S17W55 L= 31)   374     65        0 (1678/ 620-1739/ 669)  45.0 W52-W57 S19-S15
  

    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