Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1219 1835 (N53E20 L=292)   534     55        0 (1186/1815-1249/1854)  59.0 E23-E16 N50-N55
      2  644 1792 (N29E57 L=256)  1130    152        1 ( 604/1723- 699/1852) -56.3 E61-E52 N23-N35
      3  864 1814 (N37E42 L=271)   467     59        1 ( 842/1788- 895/1834)  33.7 E44-E39 N34-N39
      4  877 1725 (N31E34 L=279)   697    154        1 ( 855/1650- 908/1795) -39.8 E39-E27 N26-N36
      5  889 1441 (N13E20 L=293)  1086    286        2 ( 861/1310- 918/1580) -63.4 E25-E15 N05-N21
      6  218 1152 (S19E60 L=252)   386     88        0 ( 194/1116- 236/1200)  68.2 E62-E58 S20-S15
      7 1428  607 (S16W37 L=350)   438     60        1 (1404/ 577-1461/ 642) -63.4 W35-W39 S18-S14
      8 1529  470 (S20W52 L=  5)  1588    274        5 (1423/ 415-1633/ 550)  60.3 W47-W56 S26-S12
      9 1245  511 (S28W31 L=343)  1009    217        2 (1146/ 491-1353/ 540)  42.0 W25-W36 S31-S22
     10  514  373 (S65E24 L=289)   526     44        1 ( 489/ 357- 542/ 397)   0.0 E28-E18 S67-S62
     11  684  305 (S64W03 L=316)  1661    194        2 ( 604/ 268- 761/ 359)  10.8 E08-W14 S67-S59
     12  836  316 (S55W17 L=329)   436     54        1 ( 817/ 284- 863/ 344) -80.5 W13-W18 S58-S52
  

    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