Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  269 1452 (N23E67 L=  4)  1394    201        1 ( 237/1370- 308/1542) -71.6 E71-E64 N17-N29
      2  480 1511 (N29E46 L= 25)   671     91        2 ( 440/1469- 514/1541)  76.0 E50-E42 N26-N30
      3 1838 1293 (N23W65 L=136)   552     68        0 (1815/1253-1857/1332) -90.0 W62-W67 N20-N25
      4  281  838 (S15E53 L= 18)  2376    427       11 ( 147/ 717- 421/ 935) -19.0 E69-E40 S21-S09
      5  817  870 (S09E12 L= 59)   255    113        1 ( 772/ 841- 869/ 909)  33.7 E14-E08 S10-S05
      6  221  733 (S22E63 L=  8)  7041    469       19 ( 122/ 564- 337/ 886)  74.1 E79-E48 S34-S13
      7 1609  620 (S20W45 L=116)  5399    724       11 (1328/ 475-1828/ 860)  38.5 W25-W61 S31-S04
      8 1225  753 (S13W14 L= 85)   271     94        1 (1193/ 733-1275/ 782) -21.8 W12-W17 S14-S11
      9  765  585 (S27E15 L= 56)  1557    256        6 ( 647/ 554- 882/ 620)   6.7 E23-E06 S29-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