Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  588 1696 (N30E46 L=337)   738     99        2 ( 548/1658- 633/1727)   0.0 E50-E41 N26-N32
      2  374 1620 (N22E62 L=320)   486     69        0 ( 342/1587- 401/1646)  59.0 E63-E59 N19-N24
      3  479 1617 (N23E50 L=332)   356     67        0 ( 449/1588- 508/1641)  53.1 E52-E48 N20-N24
      4 1288 1346 (N17W11 L= 33)   321     90        0 (1260/1320-1314/1367) -26.6 W08-W12 N15-N18
      5  442  942 (S19E37 L=346)   236     47        0 ( 421/ 915- 467/ 965)  36.9 E38-E34 S20-S17
      6  974  765 (S22W01 L= 23)   949    272        0 ( 861/ 744-1113/ 792)  19.7 E05-W09 S24-S19
      7 1710  642 (S15W55 L= 78)   411     39        1 (1687/ 625-1736/ 666)  26.6 W52-W58 S16-S14
      8 1538  475 (S29W49 L= 72)   489     73        0 (1509/ 439-1571/ 508)  90.0 W47-W50 S31-S26
      9  648  351 (S57E21 L=  1)   517     33        2 ( 624/ 324- 672/ 381) -78.7 E24-E18 S59-S54
  

    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