Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  911 1785 (N60W13 L=287)   379     43        1 ( 890/1757- 930/1806)  68.2 W10-W17 N57-N62
      2 1489 1610 (N32W49 L=323)  1916    240        2 (1400/1575-1601/1645) -28.1 W41-W59 N27-N35
      3  599 1466 (N41E24 L=250)   983    207        1 ( 499/1447- 705/1490)  -9.0 E33-E14 N38-N44
      4  242 1454 (N45E70 L=204)  1045     93        2 ( 207/1423- 294/1474)   8.1 E83-E60 N42-N46
      5 1825 1267 (N03W63 L=336)   490     46        1 (1802/1243-1847/1285) -56.3 W59-W66 N01-N04
      6 1759  677 (S31W48 L=322)  2834    372        3 (1643/ 570-1869/ 837)  35.8 W36-W60 S36-S22
      7  234  625 (S07E69 L=205)  1220     98        5 ( 205/ 567- 258/ 668)  60.3 E73-E64 S11-S04
      8  646  463 (S22E37 L=237)   275     43        0 ( 624/ 439- 669/ 487) -76.0 E37-E35 S23-S19
      9 1063  275 (S44E15 L=259)  3343    490        5 ( 864/ 182-1284/ 389)  10.5 E38-W05 S48-S38
  

    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