Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  936 1610 (N39E14 L=126)  1028    175        2 ( 869/1552-1019/1654) -23.2 E20-E06 N35-N41
      2  463 1511 (N26E49 L= 91)  1654    257        6 ( 403/1403- 512/1587) -34.5 E59-E42 N19-N31
      3 1001 1390 (N24E05 L=135)   508    128        2 ( 979/1331-1022/1465)  83.7 E06-E03 N20-N29
      4  681 1377 (N20E27 L=113)   256     96        0 ( 650/1344- 724/1396)  45.0 E28-E23 N17-N21
      5  606  787 (S17E25 L=114)  2969    624        5 ( 386/ 599- 779/1000)  45.0 E46-E14 S32-S02
      6 1823  703 (S12W67 L=207)  1928    256        6 (1772/ 656-1871/ 787)  53.1 W60-W75 S14-S06
      7 1558  565 (S23W44 L=184)  2347    674        2 (1339/ 416-1697/ 708)  46.4 W30-W53 S35-S12
  

    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