Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1238 1335 (N21W20 L=195)   504    196        0 (1146/1314-1323/1358) -31.0 W14-W24 N18-N24
      2 1849  895 (S19W58 L=234)   673    118        0 (1822/ 837-1875/ 945)  45.0 W54-W62 S22-S16
      3 1709  896 (S15W43 L=218)   358     92        0 (1671/ 855-1740/ 934)  26.6 W39-W46 S16-S13
      4 1749  846 (S19W47 L=222)   385    107        1 (1708/ 820-1788/ 869) -20.6 W42-W50 S21-S17
      5 1098  744 (S11E01 L=175)   219    105        0 (1050/ 714-1139/ 765) -45.0 E03-W01 S13-S09
      6 1790  665 (S33W54 L=229)   483     42        1 (1772/ 643-1809/ 696) -45.0 W50-W57 S35-S30
      7 1106  643 (S18E02 L=173)   432    176        1 (1031/ 614-1172/ 686)  -7.1 E06-W02 S18-S15
      8 1563  571 (S33W29 L=205)  1441    226        1 (1508/ 475-1621/ 672)  53.1 W24-W34 S39-S26
      9  264  608 (S08E66 L=109)   571     70        1 ( 246/ 574- 281/ 644)  68.2 E69-E64 S10-S05
     10  617  521 (S17E37 L=139)  1044    288        0 ( 474/ 499- 733/ 548)  -5.7 E47-E27 S19-S14
     11 1028  372 (S35E14 L=161)  3134    536        8 ( 862/ 219-1194/ 543)  28.5 E34-W01 S46-S25
  

    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