Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1071 1848 (N59W45 L=163)   439     61        0 (1055/1817-1086/1879)  15.9 W39-W53 N56-N60
      2 1679 1588 (N21W69 L=187)   876     97        0 (1629/1566-1721/1605) -24.0 W63-W73 N18-N22
      3 1365 1544 (N28W38 L=156)   477     84        1 (1325/1516-1394/1573)   7.1 W33-W41 N27-N29
      4  249 1067 (N24E56 L= 62)   443     86        1 ( 230/1025- 271/1102) -68.2 E58-E53 N21-N26
      5  321 1008 (N20E49 L= 69)   672    203        2 ( 298/ 940- 383/1094)  84.3 E50-E43 N15-N25
      6 1907  993 (S19W67 L=185)   700     69        1 (1890/ 961-1925/1030)  53.1 W63-W71 S21-S17
      7 1743  914 (S18W45 L=163)   832    170        0 (1701/ 842-1794/ 977)  31.0 W40-W51 S21-S14
      8 1521  895 (S12W27 L=145)   492    112        0 (1465/ 881-1572/ 914) -15.9 W23-W30 S14-S11
      9  859  459 (S23E25 L= 93)  4472    895        8 ( 629/ 243-1116/ 707)  20.6 E57-E01 S35-S13
     10 1226  492 (S30E00 L=118)   343     97        1 (1181/ 470-1259/ 526)  21.8 E03-W02 S31-S28
  

    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