Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1317 1453 (N22W23 L=255)   682    207        0 (1221/1431-1410/1476)   4.1 W16-W30 N20-N23
      2  452 1429 (N29E37 L=195)   647    172        0 ( 393/1405- 517/1458) -20.0 E42-E31 N27-N31
      3 1013 1339 (N18W02 L=234)   229     58        0 ( 973/1324-1054/1353) -26.6 W00-W04 N16-N18
      4  294 1290 (N22E49 L=183)   328     56        0 ( 272/1262- 316/1316) -45.0 E51-E46 N19-N23
      5 1905 1273 (N07W72 L=304)   642     77        1 (1884/1254-1923/1295)  12.5 W67-W76 N05-N08
      6  263 1181 (N15E52 L=180)   335     47        2 ( 246/1152- 280/1213) -63.4 E53-E49 N13-N17
      7 1230 1137 (N04W13 L=245)   267     61        0 (1216/1096-1245/1178) -90.0 W12-W14 N01-N06
      8  281 1093 (N10E50 L=182)   889    136        1 ( 223/1054- 323/1135) -32.0 E55-E46 N07-N12
      9 1558 1069 (S03W34 L=266)   620    161        0 (1493/1013-1621/1097) -36.9 W29-W38 S06-N00
     10  186  922 (N01E62 L=170)   491     96        1 ( 159/ 899- 210/ 948)   9.5 E65-E58 N00-N01
     11 1379  882 (S13W21 L=253)   316     53        0 (1361/ 841-1395/ 915)  90.0 W19-W21 S15-S10
     12  861  862 (S10E11 L=221)   249     29        1 ( 832/ 840- 884/ 881)   0.0 E13-E09 S10-S08
     13  215  705 (S13E64 L=168)  3047    401        5 ( 147/ 613- 279/ 834) -66.8 E69-E58 S18-S04
     14 1358  720 (S23W19 L=251)   532    171        0 (1290/ 640-1432/ 800)  45.0 W14-W23 S27-S18
     15 1255  729 (S21W12 L=244)   365    126        0 (1208/ 681-1294/ 797)  45.0 W08-W14 S23-S16
     16  480  536 (S26E44 L=188)   775    182        1 ( 422/ 506- 537/ 577)  -6.3 E48-E38 S27-S23
     17 1387  388 (S47W23 L=255)  1250    207        3 (1291/ 346-1478/ 446)  14.9 W15-W31 S49-S42
  

    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