Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  981 1667 (N43W17 L=206)   549    133        1 ( 946/1625-1037/1702)  11.3 W13-W23 N40-N44
      2  437 1306 (N34E33 L=156)  1160    214        3 ( 373/1242- 530/1373)  12.1 E38-E24 N30-N36
      3  992 1294 (N19W05 L=193)   288     90        1 ( 978/1246-1009/1340)  59.0 W02-W06 N16-N21
      4  192 1250 (N36E61 L=128)   601     90        0 ( 144/1231- 237/1265) -10.6 E70-E54 N33-N36
      5 1891  990 (S21W62 L=250)  2437    295        4 (1856/ 904-1949/1118) -65.0 W57-W73 S27-S12
      6 1747  974 (S17W45 L=233)   265     64        1 (1702/ 953-1773/ 987) -26.6 W41-W47 S18-S15
      7  804  848 (S01E17 L=172)   248    123        0 ( 791/ 785- 816/ 902)  59.0 E17-E14 S04-N01
      8 1508  768 (S23W23 L=212)   331     73        0 (1488/ 728-1523/ 809) -90.0 W21-W24 S25-S19
      9  196  758 (N06E66 L=123)   510     52        2 ( 170/ 739- 226/ 784) -15.9 E69-E62 N04-N07
     10 1689  694 (S33W37 L=226)  1021    188        2 (1642/ 622-1728/ 784)  45.0 W32-W41 S37-S28
     11  826  690 (S11E19 L=169)   317     76        2 ( 787/ 663- 864/ 714) -36.9 E21-E17 S12-S08
     12  654  347 (S28E46 L=143)   437    114        0 ( 636/ 294- 672/ 399)  37.9 E49-E40 S31-S24
  

    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