Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  658 1701 (N62E07 L=245)   422     84        0 ( 634/1663- 691/1734)  35.5 E10-E01 N58-N63
      2 1275 1639 (N35W38 L=290)  1406    240        1 (1164/1581-1379/1683) -53.1 W33-W42 N28-N40
      3  441 1617 (N62E40 L=212)   435     56        1 ( 413/1597- 464/1642)  23.2 E44-E35 N59-N63
      4  317 1153 (N31E47 L=205)   817    145        1 ( 243/1128- 398/1180) -16.4 E55-E38 N27-N33
      5 1794  915 (S22W49 L=301)  1764    283        1 (1708/ 804-1883/1007)  15.3 W38-W60 S25-S18
      6 1697  724 (S30W37 L=289)   638     97        1 (1650/ 701-1750/ 746) -11.3 W32-W42 S31-S28
      7  895  432 (S25E25 L=227)   746    219        0 ( 832/ 393- 945/ 488) -33.7 E29-E21 S27-S21
      8 1545  425 (S47W21 L=273)   545    114        0 (1515/ 377-1590/ 466)  35.5 W17-W25 S51-S44
      9  561  436 (S17E47 L=205)   450    101        0 ( 508/ 416- 605/ 454) -23.2 E51-E43 S18-S14
     10  525  347 (S22E57 L=195)  2966    220        4 ( 441/ 294- 626/ 413) -45.0 E66-E50 S26-S16
     11  628  385 (S22E45 L=207)   323     60        2 ( 594/ 368- 660/ 402)   0.0 E48-E42 S23-S20
  

    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