Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1293 1398 (N25W06 L= 36)   829    250        1 (1236/1283-1364/1495) -37.6 W00-W13 N19-N29
      2  834 1310 (N07E18 L= 12)   305    141        0 ( 779/1281- 919/1336)  45.0 E20-E14 N04-N10
      3  488  907 (S24E31 L=360)   803    204        2 ( 406/ 871- 563/ 968)  38.7 E36-E25 S27-S19
      4  324  878 (S30E44 L=347)   289    101        1 ( 308/ 830- 348/ 919)  71.6 E44-E41 S32-S26
      5 1158  798 (S13W14 L= 44)   271     82        0 (1118/ 786-1208/ 812)  31.0 W11-W16 S14-S11
      6  864  730 (S25E02 L= 29)   268     86        0 ( 827/ 708- 918/ 760)   0.0 E04-W02 S25-S23
      7  540  634 (S40E23 L=  8)  1199    205        3 ( 473/ 556- 625/ 697) -19.7 E28-E14 S43-S36
      8 1450  622 (S15W37 L= 68)   519    155        0 (1417/ 552-1497/ 689) -69.4 W35-W40 S18-S10
      9  764  570 (S38E04 L= 26)   396     88        0 ( 716/ 552- 814/ 585)   8.1 E07-W00 S39-S36
     10 1286  552 (S24W30 L= 60)   469    193        0 (1198/ 538-1377/ 566)  20.0 W24-W35 S25-S20
  

    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