Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  796 1787 (N39E44 L=163)   525     74        2 ( 769/1743- 816/1820) -29.1 E47-E38 N36-N41
      2 1410 1728 (N55W10 L=217)   439     44        2 (1376/1713-1435/1748)  45.0 W06-W12 N52-N58
      3 1166 1658 (N40E09 L=198)   604     91        3 (1132/1608-1191/1704) -45.0 E12-E05 N36-N43
      4 1711 1476 (N44W44 L=251)   468     64        1 (1681/1454-1738/1502) -26.6 W39-W47 N43-N46
      5  345  921 (S24E44 L=163)   419     56        0 ( 317/ 890- 364/ 945) -26.6 E46-E41 S26-S22
      6 1535  762 (S05W38 L=245)   254     24        1 (1513/ 742-1560/ 782)  45.0 W36-W39 S06-S03
      7  410  549 (S48E38 L=169)   706     70        1 ( 379/ 513- 439/ 590)  41.2 E43-E34 S51-S44
      8  440  478 (S53E36 L=171)   773     95        2 ( 405/ 437- 474/ 513)  49.4 E42-E31 S56-S49
      9 1539  383 (S27W59 L=266)  2551    250        4 (1422/ 343-1655/ 435)  38.7 W50-W68 S32-S20
     10  527  393 (S58E26 L=181)   499     47        0 ( 502/ 368- 548/ 421)  54.5 E30-E23 S61-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