Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  373 1567 (N28E55 L=131)  1091    206        4 ( 338/1502- 416/1650) -27.6 E68-E47 N23-N34
      2 1670 1616 (N40W53 L=240)   773    114        0 (1628/1580-1715/1648)  19.3 W45-W65 N36-N43
      3  214 1357 (N14E64 L=123)  4429    412        9 ( 117/1190- 300/1480)  60.9 E75-E54 N03-N21
      4  276 1260 (N07E54 L=133)   480     89        0 ( 251/1219- 296/1315)  78.7 E55-E51 N04-N10
      5 1139  895 (S11W08 L=194)   251     99        0 (1097/ 870-1179/ 928)  36.9 W05-W09 S13-S09
      6  172  857 (S17E65 L=122)  1878    184        8 ( 133/ 783- 207/ 924)  41.6 E72-E60 S20-S12
      7  814  759 (S21E12 L=175)   232     85        0 ( 786/ 729- 850/ 795) -45.0 E13-E09 S23-S19
      8 1488  720 (S20W33 L=219)   491     96        3 (1442/ 695-1535/ 744)   0.0 W29-W36 S20-S18
      9  538  575 (S34E34 L=152)  1498    296        6 ( 448/ 471- 647/ 681)  32.5 E46-E23 S41-S27
     10 1619  530 (S30W50 L=237)   978    186        0 (1531/ 510-1699/ 559)   3.4 W40-W58 S31-S28
     11 1409  468 (S36W34 L=221)   978    181        3 (1328/ 427-1482/ 527)  41.6 W28-W38 S39-S31
  

    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