Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  884 1532 (N38E07 L=302)   546    151        0 ( 811/1493- 959/1559)  10.3 E12-E01 N35-N39
      2  556 1529 (N39E35 L=274)   396     90        1 ( 513/1504- 605/1555) -11.3 E40-E30 N37-N41
      3 1780 1421 (N22W65 L= 14)  2659    224        7 (1730/1341-1848/1519) -67.4 W60-W73 N17-N29
      4 1663 1444 (N25W52 L=  1)   341     32        0 (1646/1419-1679/1470)  76.0 W50-W54 N23-N27
      5 1470 1390 (N24W34 L=343)   238     97        0 (1425/1374-1522/1406) -23.2 W30-W37 N22-N25
      6  256 1216 (N19E57 L=252)  1138    152        2 ( 228/1146- 285/1293)  77.5 E59-E54 N14-N23
      7 1820 1263 (N12W62 L= 11)   438     50        0 (1802/1245-1838/1284) -45.0 W59-W65 N11-N13
      8  967 1150 (N13E03 L=306)   279     94        1 ( 924/1127-1016/1168)   0.0 E05-W00 N11-N13
      9  390 1055 (N09E43 L=267)   545    157        2 ( 353/1004- 452/1090)  29.7 E45-E37 N05-N11
     10  667  901 (S01E23 L=286)  1808    434        1 ( 599/ 764- 804/1076) -58.6 E26-E15 S09-N09
     11 1526  941 (S04W32 L=341)   743    207        0 (1438/ 890-1598/1012)  18.4 W25-W37 S06-N00
     12  310  753 (S10E52 L=257)   459    113        1 ( 294/ 702- 323/ 806)  80.5 E54-E50 S12-S06
     13  316  629 (S18E56 L=253)  1961    307        3 ( 247/ 545- 407/ 733) -53.1 E63-E47 S23-S11
     14  957  563 (S24E07 L=302)   507    204        0 ( 867/ 542-1056/ 593) -18.4 E12-W00 S25-S21
     15 1065  333 (S42E02 L=307)  2186    353        4 ( 969/ 251-1158/ 465)  43.2 E11-W05 S49-S32
     16  626  259 (S47E47 L=262)  1872    165        2 ( 537/ 233- 709/ 289)   7.6 E67-E35 S50-S44
  

    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