Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  622 1639 (N30E36 L=192)   306     61        0 ( 594/1623- 654/1656) -11.3 E38-E33 N28-N30
      2  767 1522 (N23E22 L=206)   307    100        1 ( 736/1495- 806/1552) -26.6 E25-E19 N21-N24
      3  380 1356 (N09E47 L=182)   323     86        0 ( 365/1312- 396/1404)  90.0 E47-E45 N06-N12
      4  454 1247 (N03E39 L=190)   353     63        0 ( 434/1207- 477/1286)  68.2 E39-E37 N00-N05
      5  761 1116 (S03E17 L=212)   185     83        0 ( 727/1089- 798/1145) -26.6 E18-E14 S03-S01
      6 1404  821 (S13W26 L=254)  1482    366        3 (1261/ 722-1542/ 979)  52.1 W18-W33 S20-S02
      7  636  778 (S24E23 L=206)   220     68        0 ( 604/ 763- 676/ 788)   0.0 E24-E19 S24-S22
      8  805  698 (S27E11 L=218)   272     91        0 ( 763/ 669- 844/ 732)  45.0 E13-E08 S29-S24
      9 1576  661 (S20W42 L=270)   338     90        0 (1545/ 616-1597/ 708)  80.5 W40-W42 S23-S17
     10  647  344 (S55E27 L=202)   459     84        3 ( 619/ 317- 680/ 370)  26.6 E30-E22 S57-S52
  

    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