Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  566 1538 (N20E42 L= 43)   356     74        0 ( 535/1521- 605/1551)  26.6 E44-E39 N19-N21
      2  414 1248 (N00E44 L= 41)   321     98        1 ( 392/1211- 444/1290) -26.6 E46-E40 S01-N02
      3 1350 1182 (N14W17 L=102)   423    144        0 (1303/1127-1408/1249) -29.1 W12-W21 N11-N17
      4 1865 1098 (N20W60 L=145)   443     92        1 (1840/1082-1892/1115)   8.1 W56-W63 N18-N21
      5  745  943 (S12E15 L= 70)   353    141        2 ( 718/ 872- 761/1009) -66.8 E17-E13 S15-S08
      6 1882  877 (N07W68 L=153)   848    146        1 (1851/ 838-1905/ 910) -16.7 W62-W73 N05-N08
      7 1517  849 (S02W34 L=119)   421     99        1 (1486/ 805-1541/ 907)  90.0 W32-W35 S04-N01
      8  977  600 (S28W06 L= 91)   853    358        2 ( 813/ 551-1159/ 642)  34.4 E03-W16 S34-S21
      9 1484  557 (S20W41 L=126)   337    110        3 (1436/ 500-1520/ 598)  76.0 W39-W42 S24-S16
     10 1617  522 (S19W53 L=138)   824    173        1 (1543/ 474-1675/ 575)  52.1 W48-W56 S23-S14
     11 1620  439 (S24W58 L=143)   443    150        0 (1582/ 418-1648/ 455)   0.0 W53-W61 S25-S22
  

    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