Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1543 1695 (N32W59 L= 26)   665     97        1 (1513/1666-1578/1714) -31.0 W55-W61 N30-N33
      2  999 1575 (N41W13 L=340)   458     79        2 ( 964/1540-1024/1614)  26.6 W09-W15 N38-N43
      3 1741 1413 (N11W60 L= 27)   412    120        2 (1716/1381-1757/1452)  68.2 W56-W63 N08-N13
      4 1014 1409 (N30W08 L=335)   283     57        2 ( 980/1392-1050/1428) -14.0 W06-W10 N28-N31
      5 1384 1408 (N21W32 L=359)   599    200        1 (1285/1394-1476/1426) -18.4 W25-W37 N17-N23
      6 1668 1369 (N11W51 L= 18)   664    144        4 (1624/1322-1714/1403) -29.7 W47-W55 N08-N13
      7 1273 1240 (N14W20 L=346)   291     36        0 (1244/1215-1300/1261) -56.3 W17-W20 N12-N15
      8  980 1176 (N17W01 L=328)   202    134        0 ( 915/1156-1044/1187) -29.7 E03-W04 N14-N18
      9 1304 1039 (N02W17 L=344)   331     89        0 (1282/ 993-1337/1087) -90.0 W15-W17 S01-N04
     10  258  851 (N10E57 L=270)   495     94        0 ( 238/ 818- 292/ 903) -63.4 E58-E53 N07-N13
     11 1415  717 (S20W17 L=344)   453    118        2 (1375/ 662-1458/ 769)  35.5 W13-W20 S22-S17
     12 1216  424 (S35E02 L=325)   347     39        0 (1191/ 401-1245/ 444) -18.4 E03-W00 S36-S32
  

    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